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CN . Abstract 

T— I ■ _ _ _ 

^ \ Flavor singlet combinations of quark operators Og = uTu + dTd + sTs contribute to many 

^ ■ important physical observables in the low energy region of QCD. Experimentally one finds 

Q ■ the values of some of these observables to be in sharp contrast to the naive (perturbative) 

^ , theoretical expectations. This indicates that non perturbative vacuum properties might 

Q^ [ play a major role in the comprehension of these phenomena. An example of such a vacuum 

G^ I contribution is the axial anomaly, which appears in the divergence of the flavor singlet axial 



c^ 



current and which is connected to the topological properties of QCD. 

From a field theoretical point of view flavor singlet matrix elements differ from non 
Oh. singlet amplitudes in the occurrence of so called disconnected insertions. These are corre- 

r-| ! lations of hadron propagators with quark- ant iquark loops or correlations between quark- 

antiquark loops, which are mediated by vacuum fluctuations. According to their respective 
flavor composition, the disconnected insertions cancel largely in non singlet processes, but 
add in flavor singlet amplitudes. 

The lattice approach provides an ideal tool to study flavor singlet phenomena. Being 
a first principle method it should be capable to uncover non perturbative vacuum contri- 
butions and to yield, on the long run, reliable results for the size of such contributions in 
QCD. 

The present article reviews the status of flavor singlet matrix element calculations in 
lattice QCD with respect to methods, results and reliability. Special emphasis is paid to 
the discussion of state of the art calculations of the pion nucleon sigma term (JjrAf j the flavor 
singlet axial coupling of the proton G\, and the rj' mass. 



X 



S. Gijsken, Flavor singlet phenomena in lattice QCD 



Contents: 



1 Introduction 3 

1.1 The Subject] 4 

1.2 The Method 6 



2 The Pion-Nucleon-Sigma Term 



10 

2.1 Operators and Correlation Functions 12 



2.3.1 Quenched QCD 19 



2.2 What has to be computed 15 



2.3 Lattice Results 19 



2.3.2 FuU QCD 23 



2.4 Discussion 28 



3 The Flavor Singlet Axial Coupling of the Proton 



3.1 Preparation 



30 

31 

3.2 Direct Method) ■ ■ ■ ■ 35 

3.2.1 RenormalizatioE| 37 

3.3 Topological Method | 39 

43 



3.4 Lattice Results 



3.4.1 Quenched QCD 43 



3.4.2 FuU QCD 47 



3.5 Summary 51 



4 The Mass of the ?/ Meson and the U(l) Problem 



52 

4.1 The Problem! 52 



4.2 Lattice Technique 54 

O Lattice Results! 57 



4.4 Outlook 63 



5 Final Remarks 



65 



S. Giisken, Flavor singlet phenomena in lattice QCD 3 

1 Introduction 

Quantum field theories are today by far tlie most successful candidates in describing the 
fundamental forces of nature, and there is still no evidence that the standard model - a 
Lorentz and gauge invariant quantum field theory with gauge group SU{3)c x SU{2)w x 
U{1)y - fails in the description of electroweak and strong phenomena. 

One of the most fascinating consequences of such theories is that the physical vacuum 
is not vacuous. It rather behaves like a Lorentz invariant medium, which is made up from 
ground state fluctuations of quantum fields. This property, which is absent in classical field 
theories and also in quantum theory, implies a tremendous spectrum of new and partly 
curious physics phenomena, whose range has still not been fully explored yet. 

A particle brought into such a vacuum, interacts with it. Excitation, polarization 
and anti-polarization effects occur and give rise to quite singular phenomena such as the 
decay of particles into classically forbidden channels, e.g. the decay $ — > Svr, the running 
of the classically constant coupling with the distance of two particles or the property of 
confinement of strongly interacting particles |]l|. 

Due to its internal dynamics, the vacuum itself may exhibit a non-trivial structure. 
Symmetries, which are present on the classical level, can be destroyed by quantum fluctu- 
ations. 

A qualitatively well understood consequence of such a spontaneous symmetry breaking 
is the appearance of massless (Nambu-Goldstone) bosons in case of a broken global symme- 
try P,||, and of (massive) gauge bosons 0-0] in case of a broken (local) gauge symmetry. 
The former mechanism, applied to the chiral SU{3) symmetry, explains the smallness of 
the masses of the pseudo scalar meson octet members. The latter allows to include massive 
W and Z gauge bosons into the theory of electroweak interactions. 

A third and much more hidden way of symmetry breaking occurs, if the quantum 
fluctuations form an anomaly, i.e. if the conservation of a classical current gets spoiled by 
vacuum effects. In the standard model, this is known to happen to the axial vector current 
P, ^]. On the one hand the implications of such symmetry breaking can be disastrous, 
if the anomalous current couples to a gauge field. In that case the renormalizability of 
the theory, and thus its physics content, gets lost |T^ unless the gauge properties and the 
particle content of the theory are arranged such that the anomalous contributions cancel 
each other. Fortunately, the standard model is constructed in this way. On the other 
hand, if the anomalous current does not couple to a gauge field, this 'dynamical' symmetry 
breaking supplies a new handle to understand quite a number of physics phenomena. 
Besides the successful calculation of the decay ttq -^ 77, 0, it has been shown ||lT| that 
the anomaly of the flavor singlet axial vector current is related to a non-trivial topological 
structure of the vacuum, which in turn can be used to explain qualitatively the large mass 
of the 7]' meson |l^,|l3l. Moreover, the unexpectedly small value of the flavor singlet axial 



vector coupling of the proton G\, which lead to the so called proton spin crisis, is also 
closely connected to the influence of this anomaly [p!4| -[T^ . 

We do not attempt to discuss all the implications of the quantum vacuum here. We 
rather want to concentrate on a certain class of phenomena in the deeply non-perturbative 
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Figure 1: Connected (a) and disconnected (h) insertions of a proton interacting with a 
current j. Please note that all quark lines, including the quark loop, are connected by 
infinitely many gluon lines and virtual quark loops. 

regime of the quantum field theory of strong interactions (QCD), which are expected to 
be determined to a large extend by vacuum properties: The flavor singlet phenomena. 



1.1 The Subject 

Suppose a (composite) particle, say a proton, interacts with a current j. Then quantum 
field theory tells us that, in addition to the direct coupling, a second process contributes, 
which is the interaction of the current j with a vacuum quark-antiquark loop in the field 
of the proton. Both contributions, the connected (direct) and the disconnected^ (vacuum) 
insertion, are shown schematically in fig.|l]. 

One would expect however, that the disconnected insertion is small compared to the 
connected one if the current j has a flavor non-singlet structure. Consider for example the 
flavor triplet and flavor octet currents 



J3 = uTu — dVd , jg = uTu + dVd — 2sTs , 
where F denotes a combination of Dirac matrices, and the electromagnetic current 



Jem 



2/3m7^m - l/M-i^d - l/3s7''s . 



(1) 



(2) 



Assuming an approximate flavor symmetry of the vacuum, the contributions from each 
quark flavor u, d, s to the disconnected insertion would largely cancel each other. 
The situation changes drastically if we choose j to be a flavor singlet current 



ji = uVu + dVd + sVs . 



(3) 



In this case, the single quark contributions simply add up, and there is no a priori reason 
why the disconnected part should be small any more. Thus, flavor singlet processes are 
promising candidates to study vacuum polarization and vacuum structure effects. 

^Note that these 'disconnected' insertions are stiU connected in the field theoretical sense. There are 
however no valence quark lines which connect the quark-loop and the Proton propagator. 
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So far we have introduced flavor singlet phenomena solely as generic constructs. Of 
course, to be able to extract physics information from such processes one has to know to 
which physical situations they are related, whether they can be characterized by observ- 
ables which are accessible to experiment and, last not least, how the diagrams in fig.|l] can 
be calculated reliably in QCD. Only then, by comparison of the theoretical expectation 
with the experimental measurement, flavor singlet phenomena can help to scrutinize the 
quantum field theoretical picture of the vacuum. 

The question of how and to what precision flavor singlet matrix elements can be calcu- 
lated is one of the major issues of this review. As mentioned above, we will concentrate on 
flavor singlet phenomena in the deeply non-perturbative regime of QCD, i.e. at momentum 
transfers below 0{Aqcd)- The only known technique, which allows to solve QCD in this 
regime without additional assumptions is the method of lattice QCD. Thus we will use this 
method here. 

Flavor singlet processes contribute to many physical observables, the most important 
ones being related to the chiral structure of the strong interaction. A scalar current 
(ji(r = 1)), multiplied by a (quark) mass has the structure of a mass term in the QCD 
Lagrange density. Indeed, it turns out that ji(r = 1), when inserted between hadron 
ground states, just measures the explicit breaking of the chiral symmetry in QCD. The 
quantity mq{N\uu + dd\N) is known as the pion-nucleon-sigma term (Tj^n- Its value can 



be extracted from vr A^ scattering data |T^ , although the analysis is rather indirect and 
requires much theoretical input [0. 

A very interesting situation arises if we choose ji to be a axial vector current (F = 7^75) . 
The conservation of this current is spoiled by the axial anomaly, and one would expect 
to track the effect of this anomaly in the disconnected insertion of fig.0. In fact, the 
corresponding amplitude, connected + disconnected, just defines the axial flavor singlet 
coupling of the proton G\, which is extracted from the flrst moment of the spin dependent 
proton structure function g^ . Naively, G\ can be interpreted as the fraction of the the 
proton spin carried by quarks. It was flrst measured from polarized muon proton scattering 



by EMC [|19|. The result, which has been consolidated by several succeeding experiments 
pO|^4[ revealed a much smaller singlet coupling than expected [^]. One might wonder 
of course whether this result can be derived directly from QCD by a calculation of the 
diagrams in flg.[l|. 

The flavor singlet tensor charge of a nucleon, Qt = (A^|ji(F = <J^ul5)\N) can be 
extracted from the flrst moment of the chirality violating structure function hi of a nucleon. 
Since the chiral anomaly is absent here, one would not expect an anomalous result as found 
for G\. Although experimental measurements of Qt are not available yet, it is interesting 
to investigate by a QCD calculation whether this expectation is justifled. 

Closely related to the axial vector current is the pseudo scalar current (F = 75), as 
the divergence of the former is connected to the latter by the axial Ward identity. Thus 
we expect to perceive the influence of the axial anomaly also if we choose ji to be a 
pseudo scalar current. The prominent application of this choice is the U{1) problem of 
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Figure 2: Connected (a) and disconnected (h) contributions to the tj' propagator. 

QCD, i.e. large mass of the 'would be ' Goldstone boson rj'. We mentioned above that 
a deep qualitative understanding in terms of the axial anomaly in conjunction with a 
topological structure of the vacuum has been achieved already, but of course one would 
like to calculate the r]' mass directly from QCD. The graphs which have to be evaluated 
for this purpose are slightly different from those given in fig.|l|. As can be seen from fig.^, 
the connected contribution is simply the standard propagator of a pseudo scalar particle. 
The disconnected part, which contributes in case of the flavor singlet combination, consists 
of the correlation of two (pseudo scalar) vacuum loops. Clearly, one expects that the 
contribution of this correlation gives rise to the large rj' mass. 

To summerize, flavor singlet phenomena are closely connected with the quantum struc- 
ture of the vacuum. Thus, a strong motivation to study these processes is to learn about 
this structure, and to contribute to our knowledge about whether and how the standard 
model concept of the quantum vacuum is realized in nature. 

In the following chapters we will discuss the application of lattice QCD to the deter- 
mination of (T^AT, G\, and to the U{1) problem^. 

1.2 The Method 

Although there has been tremendous progress in the qualitative understanding of vacuum 
properties over the last 30 years, a precise quantitative matching of experimental results 
with theoretical calculations is still missing in many cases. The reason is mainly not a 
lack of experimental results, but rather lies in the fact that most of the QCD vacuum 
phenomena cannot be treated within the framework of standard perturbation theory. The 
only method known today which is capable to treat non-perturbative phenomena from 
first principles, i.e. from the definition of the theory itself, is the method of lattice gauge 
theory ^^. We describe here only some important properties of lattice gauge theory which 



will enable the reader to follow the line of presentation of this review, and to develop a 
feeling for the reliability of the results. 

Lattice QCD was invented by K. Wilson [^ about 3 decades ago. The basic inten- 
tion was to setup a gauge invariant and non perturbative regularization scheme for QCD. 
This can be achieved if one puts the path functional [^, which in principle contains all 



^The results of a pioneering lattice calculation of the tensor charge Qt can be found in ref. [p6[ 
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information of a given quantum field theory, on a space-time grid with a lattice constant 
a and a lattice extension L x a. The reduction from continuous to discrete space-time 
instantly introduces an ultraviolet cutoff, which removes possible singularities due to high 
momentum limits. Since the action density can be formulated gauge invariantly on that 
grid, this procedure defines a gauge invariant regularization scheme. 

The second important ingredient, which makes quantum field theory look like statistical 
mechanics and which opens the door to a numerical treatment is the Wick rotation, i.e. 
the change from Minkowski metric of the space-time grid to Euclidian metric by a rotation 
of 90 degrees of the time axis in the complex plane. Fortunately all (time independent) 
properties of the theory remain untouched by this manipulation. The great advantage is, 
that the time dependence of n-point correlation functions is changed from an oscillatory 
to an exponential behavior. This enables for the numerical evaluation of the path integral. 

With these settings, n-point functions of (arbitrary) operators Oi, O2, ■ ■ ■ On can be 
expressed by 

{Oi{xi)02{x2)...On{Xn))= (4) 

Z-' J d[U]d[ilj]d[ilj] {Oi(x,)02(x2) . . . o„(x„)e-s^(^''^'^"'^'^«)-s-(^'«')} , 

where 

Z-^ = f d[U]d[il;]d[ilj]e-^f-^' . (5) 

The fermion fields ip and ip are anti-commuting Grassmann variables. The gluonic de- 
grees of freedom are represented by the gauge-links U^{x). Those are SU{^) color ma- 
trices, defined as the parallel transporters, which gauge-invariantly connect the fermion 
fields at lattice point x with the fermion fields at the neighboring point in the direction 
11. Sf{U,ip,'ijj, M,gQ) and Sg{U,go) denote the fermionic and gluonic parts of the (lattice 
discretized) QCD action. 

Note that there are only two free parameters in eq.(^, the bare coupling go and the 
bare quark mass^M. All dimensionful quantities are given in units of the lattice spacing a. 
We emphasize that a is not a free parameter: At a given choice of f^o? M can be tuned such 
that the lattice value of a dimensionless quantity, say the ratio M.^/Mp, corresponds to its 
physical value. The size of a is then determined by the matching condition^ Mp = nip a. 

The first step in the evaluation of the path integral can be done analytically. Since Sf 
has the generic form 

Sf = ^M{U)ij , (6) 



•^For clarity we consider here only one (light) quark mass. Although the extension to more (non- 
degenerate) quark masses is in principle straightforward, one faces severe problems in the implementation 
of the Monte Carlo process. Most (unquenched) lattice calculations are therefore done with (two or four) 
degenerate quarks. 

^Lattice values are indicated by capital letters, physical values by small letters. 
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the fermionic degrees of freedom can be integrated out |Q . Eq. (|^) then reads 

(Oi(x,)02(x2) . . . 0„(x„)) = Z-' I d[U] det(M)F([7, M-^([7))e-^- . (7) 

F is a known function, which depends on the choice of the operators O,. It can be expressed 
in terms of the inverse fermionic matrix M~^ and the gauge hnks U. For example, for the 
quark propagator, Oi{xi)02ix2) = V^(xj, a, a)V'(x2, a', a'), F is given by M~^\f^.^^,^,p,, 
where b, b', /3, /?' denote color and Dirac indices. 

To proceed further in the evaluation of eq.(^ one relies on numerical importance sam- 
pling methods @,|35|, which generate 'gauge configurations' [U] with a probability 



P[U] oc det{M[U]) X exp{-Sg[U]) . (8) 

The expectation value {Oi{xi)02{x2) . . . On(xn)) then is just the average of F[U] over all 
gauge configurations [U]. Note that all results achieved with such a Monte Carlo proce- 
dure pick up a statistical uncertainty, since one is always limited to a finite sample of gauge 
configurations. Of course, this uncertainty can be reduced systematically by increasing the 
size of the statistical sample. 

The most computer time consuming part of this stochastic approach is the generation 
of a statistically representative sample of configurations. Especially the calculation of 
det(M[f/]), which is a highly non local quantity, slows down the Monte Carlo process 
considerably, and it is only in these days that computers become powerful enough to 
create statistically significant ensembles of [U] on acceptably large lattices. Therefore, in 
the past, the bulk of lattice results has been achieved in the quenched approximation, i.e. 
by setting the fermion determinant to the fixed value det(M) = 1. In terms of physics, 
'quenching' corresponds to a neglect of fermionic vacuum loops in the process of generating 
gauge configurations. 

Naively one might expect that most of the quenching effects can be accounted for by 
a redefinition of the coupling qq. However, this expectation could be misleading when 
one studies flavor singlet phenomena. As explained above, fermionic vacuum loops play a 
crucial role in this case. Allowing for such loops in the construction of n point functions, but 
neglecting them in the creation of vacuum configurations could induce a systematic bias. 
The size of this bias can only be estimated by a comparison of quenched and unquenched 
results. This will be one important topic of this review. 

We did not yet discuss the explicit forms of Sf and Sg. Clearly there is a lot a freedom 
in the definition of these quantities, as they have to coincide with the QCD action only 
in the continuum limit, i.e. at a — »• and L — > oo. Of course one tries to construct 
the discretized action density such that as many continuum symmetries as possible are 
preserved and discretization errors are kept as small as possible. However, due to the 
so called fermion doubling problem p0|^2|, this turned out to be a highly non-trivial 
task. Two different prescriptions, one proposed by Wilson PS|,PB| and one by Kogut and 
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Susskind |^^ , have been mainly used for Sj over the recent years. With the Wilson action, 



discretization errors occur, together with an explicit breaking of the chiral symmetry, 
at 0{a) in the lattice spacing. The Kogut-Susskind method preserves the U{1) chiral 
symmetry, but mixes internal (spin) and external (space) degrees of freedom. In that 
scheme discretization errors contribute at 0{a^). 

Much effort has been invested to reduce the discretization errors of the Wilson ac- 
tion from 0{a) to 0{a'^). The most promising candidates to achieve this goal are by 



now the semi-perturbatively |39| and the non-perturbatively |3H] improved versions of the 



Sheikholeslami-Wohlert action ||40|| . 

The lattice definition of Sg is much less problematic. The Wilson prescription 
which exhibits discretization errors at 0(0^) is still mostly used, but improved versions [^ 
{0{a^)) are also available. 

In order to make contact with continuum physics, the parameters Qq and M have to be 
tuned such that the cutoff a~^ is removed while keeping the quark-mass m at its physical 
value. In QCD this means that, because of asymptotic freedom, one has to take the limits 
qq —>■ 0, M = ma —* 0. Of course, to keep the physical size La of the lattice finite, one 
has to increase L oc 1/a in this procedure. 

In practice, one takes these limits by an extrapolation of the results from several sets 
of parameters {qq, M, L) to the continuum. 

The question of how close to the continuum limit these sets of parameters can be 
chosen, is, after all, a question of computer power. In terms of statistical physics, the 
continuum limit is located at a phase transition point of second (or higher) order. Close to 
this point all correlation lengths of the system start to grow exponentially and the Monte 
Carlo process of generating decorrelated configurations slows down more and more. 

Current resources have enabled for Monte Carlo simulations with lattice specifications 
a ^ 0.05 fm , L/a ^ 3 fm , m ^ l/2mstrange in the quenched case, and a ^ 0.1 fm , L/a ~ 
2fm,m^ rristrange for full QCD Q. 



The question of how close to the continuum limit the parameter sets have to be chosen to 
allow for a reliable extrapolation, is difficult to answer in general. Although the functional 
dependence of the lattice results on M and a is well known close to their respective limits[| it 
is not a priori clear, at which values of M and a the asymptotic behavior sets in. Therefore 
one has to check the onset of asymptotics carefully with the lattice data in each case. 
Clearly, the reliability of the continuum extrapolation improves successively as more and 
more lattice results become available. 



^ Close to the continuum limit the M dependence can then be calculated by chiral perturbation theory 
51 and the a dependence is proportional to the order in a at which discretization errors occur. 
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2 The Pion-Nucleon-Sigma Term 

The Pion-Nucleon-Sigma term a^^N is defined as the double commutator of the axial charge 
with the Hamilton density, sandwiched between nucleon states at momentum zero: 

cTnN = ^(iV|[Q'(0), [Q\0),nm\N) , (9) 

where M is the nucleon mass and 7i is the QCD Hamilton density. The axial charge is 
given by 

A^{xo, x)d^x , with J^ = q'yo—q , and q = {u, d, s) . (10) 

Since ct^at contains the time derivate Q^ = [Q^,7{\ of the axial charge, it supplies a direct 
measure of the explicit chiral symmetry breaking of QCD. If the chiral symmetry would 
be conserved, the axial charge would be also conserved {Q^ = 0) and ct^at would vanish. 

In order to exploit eq.(P) one uses the assumption that the chiral symmetry breaking 
term of the Hamilton density is given by ^, ^ 



7i' = rriuUU + niddd + niaSS . (11) 

Neglecting, in addition, SU(2) flavor symmetry breaking effects, i.e. assuming ruu = rrid, 
eq.(H) takes the form 

cr.N = ^"^^^^^{Nluu + dd\N) . (12) 

Note that the matrix element on the r.h.s. exhibits a {SU(2)) flavor singlet structure, as 
discussed in the previous section. 

The determination of the physical value of aT^N from vrA^ scattering experiments requires 
a lot of theoretical input. The isospin even wN scattering amplitude D~^ is related to 
a^N at the (unphysical) Cheng-Dashen pointQ by 



S = F^D+{2ml) ,J:^a^N for m,, m^ ^ . (13) 

Since a^^N is deflned at pion momentum transfer t = 0, a number of successive extrap- 
olations is necessary to extract its value from the experimental data: One flrst has to 



extrapolate D~^ from the physical region t < to the Cheng-Dashen point [^. The most 
recent analysis of experimental data leads to a.^i\i{2m'^) ~ 60MeV |jTB[, where a (small) 
correction term A/j = S — ajrAr(2m^) due to the non vanishing quark masses has been 
included. 

The second step is to calculate the shift A^- = o"7rAr(2m^) — a-^NiO) of cr,rAf when going 
from the Cheng-Dashen point to t = 0. Gasser, Leutwyler and Sainio [l^ have performed 



^D^ is the scattering amplitude with the born term removed. 

^The Cheng-Dashen point is defined by q^ = g'^ = 277i^, t = 2to^, z^ = i^b = 0. q,q' are the momenta 

of the incoming and outgoing pion, t, i^, i^t are given by i = {q^ — g^^, i' = -jf, vb — ~~^w^- Qt^ ^^'^ 
P^ denote the pion and nucleon momentum sums: Q^ = 5 ('Zm ^ ?«) ^^^'^ ^p ^ 2 (^'m ^ Pu)- 
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a careful analysis of this shift by means of dispersive relations and chiral perturbation 
theory. They find Ao- — 15MeV, and, correspondingly 

a^N ^ 45MeV . (14) 

A similar analysis has been done in the framework of heavy baryon chiral perturbation the- 
ory by the authors of ref. [BDl. They find aT^N = 48±10MeV, which is consistent with eq.lTl. 



Additional information about u^at can be extracted from the baryon octet mass split- 
tings. On the one hand, the octet quantity 

ao = ^^^^^^^{N\uu + dd-2ss\N) (15) 

is related to a-„N by 

ao 2{N\ss\N) 

""-^-T^y^y- {N\uu + dd\N) ■ ^^^^ 

Thus, under the assumption that the strange quark contribution to the nucleon mass is 
small, (A^|ss|iV) ~ 0, one would expect a.„N — o'o- 

On the other hand, ctq is directly related to the nucleon mass shift caused by explicit 
SU{'i) flavor symmetry breaking. To first order in m^ — l/2(m„ + m^) one finds 

/^MM = -{^^^^^^^^^^^-m)i{N\uu + dd-2ss\N) , (17) 

or, for the baryon octet mass splittings in general 

AMb = {B\c8Us\B) = atriBusB) + ptr{BBus) . (18) 

B on the r.h.s. denotes the baryon matrix, and the octet quantities Cg and Ug are defined 

as 

1 / m„ + Hid \ 1 



Cg = ^ f mA , Us = —i={uu + dd - 2ss) . (19) 

From the analysis of the mass splittings one naively extracts o"o — 25MeV. It has been 
pointed out however |J5l]|, that the corrections due to terms oc (m^ — l/2(m„ + m^))^ 
enhance this value on a (20 — 30)% level. Including these terms one finds 



ao = 35 ± 5MeV . (20) 



Comparing eq.(|lj) with these findings one concludes that 

y- 0.2 -0.4. (21) 

Let us try to interpret this result in terms of connected and disconnected insertions. 
The term (A^|ss|P) receives a contribution only from the disconnected part of fig. [^, since 
the nucleon does not contain a valence strange quark. Thus, the result y 7^ is direct 
evidence, that disconnected diagrams really contribute to a-^N- 
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To get an idea of what eq.(^) can tell about connected and disconnected insertions, 



we reformulate y in terms of the ratio of disconnected to connected contributions 

^ {N\uu + dd\N)d,,, 
"" {N\uu + dd\N),on ^ ' 

and of the flavor symmetry breaking parameter a^fe of the disconnected insertions 

(JVI^^IAT) 

where qi denotes a light quark. This yields 

Thus, naively assuming flavor symmetry of the disconnected (vacuum) parts [agb = 1), 
one would need Rdc — 0.5 to achieve y ~ 0.3, i.e. the connected contribution to a^N 
should dominate by a factor of 2. On the other hand, if one would assume connected 
and disconnected parts to contribute equally to (T^at, one would need a strong symmetry 
breaking, asb — 0.6, to explain the experimental estimate of eg. (pi]). 

Of course, these assumptions are speculative and an accurate determination of the 
matrix elements in eq's.(^,^) is highly needed. It will be interesting to see, to what 
extent state of the art lattice calculations can contribute to this task. 

2.1 Operators and Correlation Functions 

We are interested in a computation of matrix elements (A^|gg|A^). The path integral ap- 
proach, c.f. eq. (H), tells us how to determine correlation functions of operators. In this 
paragraph we explain which correlation functions are needed and how they can be be com- 
bined to obtain the required information. 

Consider the correlation functions 

Gnn{xo) = E(^(0' 0)^^(^' ^o)) (25) 

X 

and 

GnM^o) = EE(^^(0.0)^(^'?/o)g(y,Z/o)iV(x,a;o)) . (26) 

^ y,yo 

Here, N denotes the creation operator of a nucleon-like particle, i.e. a particle with the 
nucleon quantum numbers, which is however not necessarily in its ground state. Such an 
operator is most easily expressed in terms of the fermion fields^ u and d p3[p| 

iV^ = EKC'75rf')u;e'^''''^ (27) 

a,b.c 



^We write down here the expression for the proton. For the neutron, u and d have to be interchanged. 
^For the time being expHcit forms of operators and correlation functions in terms of fermion fields and 



fermion propagators will be valid only for the Wilson quark action. 
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C is the charge conjugation operator, Greek letters denote Dirac indices, and small Latin 
letters are used for color indices. The time (xq) dependence of the momentum zero nucleon 
propagator Gnn is well known. It has the formQ 

Gnn{xo) = E l(0|^V)r e""^'""" "^^^ l(0|A^W)r e"*''^"" . (28) 

n 

Thus, although Gnn contains all (radial) excitations of the nucleon, the ground state is 
isolated in the limit of large xq. 

We now turn to the time dependence of the 'insertion propagator' GjsfqN- Since the 
term m^gg is contained is the QCD action, Gnqn can be constructed by deriving Gnn 
with respect to the quark mass. Because of 



{N\0,0)N{x,Xo)) = Z-^ / d[U]d[ij]d[ij]N^{0,0)N{x,Xo)e-^f^'^^^'''^''^^ (29) 

it follows that 

{N\0,0)J2{qq){y,yo)N{x,Xo))c 

Gn,n{xo) - {N^{6,0)N{x,Xo)){J2{qqMyo)) 

y,yo 

Z''^ f d[U]d[tlj]d[tlj]N^{6, 0)N{x, Xo)e-^f~^^ 



y,vo 



dlTiq 



drUq 



^ (iVt(0,0)iV(x,xo)). (30) 



Using eq. (pHf ) one finally gets in the limit of large Xo0 



jSC/M/ N _ ^AfgAf(^o) /V^/-N/-„ N\ A , ^^^ 



rS^^\x,) ^ ^l^Ji^ipi - (Y^iqqM y,)) =A + ^Xo , (31) 



Gnn{xo) T-; ' dm. 



y,yo 



1 



with A = —d/dmq{\n\{0\N^N)\'^). Note that it is the connected correlation function, in 
the field theoretical meaning 

(iV(gg)iVt), = (iV(gg)iVt) - (NN^) (qq) , (32) 

which contributes to eq. (^Tj). Thus, the lattice value of chiral condensate (qq) is explicitely 
subtracted^. 

"'^^For the sake of clarity we will always assume a lattice with infinitely large extension in time. The 
corrections due to finite time extension are well known, but they will not be discussed here. 

^^This summation method has been introduced by Maiani et al. |54|. 

'^^As the latter diverges cubicly with the cutoff a~^ Isl], this could lead to numerical problems for lattice 
simulations close to the continuum limit. 
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The next step is to determine the derivative dMN/drUq in terms of the matrix element 



{N\qq\N). The Feynman-Hellmann theorem ||56| states that 



— - (^) (33) 

for any parameter A on which the Hamiltonian Ti depends. Thus, since the exphcitly quark 
mass dependent part of 7i is given by rrigqq one finds 

^ = wmm . (34) 

One procedure to determine the relevant matrix element is therefore to calculate the corre- 
lation functions G^jy and Gjyg^ by solving the corresponding (lattice) path integrals, and 
then to extract the slope of the ratio R^^^ . We call this technique summation m,ethod. 

Clearly, according to eq. (@), one can determine {N\qq\N) also by a direct lattice 
calculation of Mi\i{mq) at several quark masses mg. However, this would require to derive 
Mj\i{mg) numerically with respect to rUg, which in practice leads to a significant loss of 
statistical accuracy. We will come back to this point later on. Note that this method does 
not discriminate between connected and disconnected insertions, and it cannot be used to 
calculate the y parameter, eq. (|1^). Applied in the quenched approximation, it yields only 
the connected part of a-^N- 

One could argue, that the requirement of the nucleon being in its ground state might not 
be strictly fulfilled for the 3-point correlator Gnqn- Although xq can be made large, there 
are always short (time) distance contributions from the correlation of the scalar insertion 
{qq){yo) with the proton, since G^qN contains the sum over i/q. To check on the influence 
of such contributions we sketch here briefly a third method to compute {N\qq\N). This 
method, which we call plateau method, uses the correlation function 

CNqN{xo, yo) = J2J2(^^(^^ o)^(^' yo)qiy, yo)N{x, xo)) , (35) 

X y 

in which we do not sum over y^. The connection to the matrix element (iV|gg|A^) can 
be most easily obtained in the (Euclidian) continuum]^. The prescriptions for shifts of 
operators O in space and time are given by 

C(0, xo) = e^^^O{x, xo)e-*^^" and C»(x, 0) = e-^^oC»(x, Xo)e^^o . (36) 

V and Ti denote the momentum and Hamilton operators. Inserting complete sets of energy 
eigenstates \n) 

Y^\n){n\ = l (37) 



^■^The corresponding calculation on the lattice would have to be performed in the somewhat cumbersome 
transfer matrix formalism [p^ . The result is however identical if we assume an infinitely extended lattice 
in time. 
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to the right and to the left of q{y,yo)q{y,yo) in Cnqn, shifting all operators to (0,0), and 
carrying out the spatial sums one finally gets for the (field theoretically connected) part 

CNgN{xo, yo) - Gnn{Y1 ^liV' yo)) 

y 

= J2{'^\N^\n){n\uu\n'){n'\N\0)e-^-y°e-^"'^''"-y°^ . (38) 

•n,n' 

The ground state \n) = \n') = \N) is isolated only if both time separations, (xq — yo) and 
yo, become large. 

Note that eq. ( ^8]) is independent of yo in the large (time) distance limit. In this region, 
the ratio 

i?^^^(xo, yo) = ^ - (E 9^(^> yo)) = (mmlN) (39) 

is just the matrix element of interest. Thus, instead of looking at a linear slope, one has 
to identify a plateau region in (xo,i/o)- 

The plateau method is clearly superior to the summation method with respect to a reli- 
able ground state identification. However, the ratio R^^^ might be more noisy than R^^^ , 
as one refrains from accumulating statistics by summing over yo- We will see below, that 
statistical fluctuations are indeed a major problem in the determination of disconnected 
contributions. In this case, even the summation method might fail, since contributions from 
the region yo ^ Xo occur, which do not add to the signal but largely enhance the noise. 

The plateau accumulation method (PAM) tries to avoid the drawbacks of the methods 
described above. It is defined by 

R^'^^'ixo, Axo, Ax',) = E ° i?^^^(xo, yo) , (40) 

yo=Axo 

with 1 < Axq, Axq < Xo- The asymptotic time dependence is given by 

i?^^*^(xo, Axo, Ax'q) = A+ {N\qq\N){xo - Ax', - Aa;o) • (41) 

With PAM the summation is restricted to the region in y, where ground state behavior of 
the 3-point correlator can be expected. Clearly, to ensure ground state behavior, one has 
to vary Axq, Axq, and to compare the results. 

In the next paragraph we will briefly explain how the correlations functions Gat at, Gnqn 
and CjyqN can be determined on the lattice. 

2.2 What has to be computed 

Assume that we would have been able to create a statistically significant set of gauge con- 
figurations [U] by some fancy Monte Carlo algorithm. Then, to evaluate the path integral 
according to eq. (|^, one 'simply' determines the function F{U,M~^{U)) for each of the 
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correlation functions Gnn, Gnqn and Cnqn, calculates its values on each gauge configu- 
ration, and finally averages over all configurations. The purpose of this paragraph is to 
construct F for the above correlations functions, and to explain how it can be calculated 
numerically on each gauge configuration. 

The key quantity, from which all our correlation functions can be constructed, is the 
propagator A, of a quark q, created at (x, xq) and annihilated at {y, yo) 



Ag(f , xo; y, yo) = g(x, xo)q{y, yo) 



(42) 



Once we know the corresponding function F/^ , it is just a matter of combinatorics to 
determine Fg^^^, Fgj^^j^ and Fcj^^^^. 

F^^ can be most easily obtained in the generating functional formalism |27,33]. One 
finds 

Fa,{x, Xo; y, yo) = M~^{y, yo] f, xq) , (43) 

where M is the fermion matrix, c.f. eq. (^. If we insert the nucleon operator, eq. (p^), into 
26|), and (|35|) and contract q and q in all possible ways, we can express Fq^^, 



eq's. (|2 
variables 



Fc^qj^ and Fcj^^j^ in terms of M ^. Note that q and q are anti-commuting Grassmann 



{q{x),q{y)} = {q{x),q{y)} = {q{x),q{y)} = 



(44) 



We do not give here the explicit expressions for all these functions as those are quite 
lengthy. We merely want to discuss their structure in order to reveal how the different 
parts can be computed. 

The propagator of a nucleon that travels from the source (0, 0) to the sink {x, Xq) 
contains the product of three quark propagators 



G'g,{0;x) = {F^':l{0;x)) 

= ((M-i)^'^'(x; 0) [C^,M^\x; 0)^75]"'"' {M'')^'^' {x; 0) 
+ (M-i)^'"'(x;0) [C75M,-i(x; 0)^75]°'"' (M„-i)"'^'(x; 0)) . 



(45) 



Here we have suppressed the color structure on the r.h.s.. With the source point (0,0) 
being fixed, the computational problem in exploiting this equation is just to determine one 
rowQ of the inverse of the fermion matrix, i.e. to solve 



M,.,,FaMx')=6. 



x,0 



(46) 



Although M is a huge matrix of size [N^ x Nt x 3 x 4 x 2)^, the solution of this equation 
is not a severe problem^ if one uses fast iterative inverters |]57|-^|. 



^''in fact, as each fermion has 3 color and 4 Dirac degrees of freedom, one has two compute 12 rows of 
^^This is true for typical lattice sizes and quark masses used in current simulations. 
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We now turn to the calculation of Gnqn- The application of the contraction procedure 
reveals two important differences to eq. (^51) . First of all, there is an additional term 
FGNN{,0]x)J2yTr{M~^{y;y)). This term can be identified immediately with the discon- 
nected contribution to Gnqn, since M~^{y,y) describes just the quark propagation from 
y ^ y. Thus, 

G%l% = {FG^,{0;x)Y.TriM-\y;y)). (47) 

y 

In our graphical language, G'}^^% is represented by fig. |l](b). The computational expense 
to calculate this contribution is huge. We will come back to this point below. 

The second difference concerns the connected part of Gnqn- Since we have inserted the 
operator J2y{(l(l){y), one of the nucleon valence quarks, started at the source 0, will always 
meet the 'inserted' quark at y, where the (scalar) interaction is located, and then travel to 
its destination at x. Thus, the propagator of this quark is modified 

M-\x,0) -. M-\x,0) = Y.M;\x,y)M;\y,0) . (48) 

y 

Graphically, such a modified quark propagator is represented by the upper quark line of 
fig.S(a). 

Therefore, GnqN can be written in the form of eq. (|i5|), but with one quark propagator 
replaced by a modified propagator. Choosing the notation Gnn = ^{M~^ ^ Mj^ ^ M~^) 
one gets in the case of an w-quark insertion 

GT,n{u) = HM-\ M^\ M-i) + T{M-\ M^\ M-') , (49) 

and for the d-quark insertion 

GT.Nid) = nM-\ M,\ M-') . (50) 

Clearly, there is no such contribution for the s-quark. 

The numerical computation of M~^ can be done in a two step procedure ||5^. One first 
determines M~^{x,0) according to eq. (1^). Then, because of 



M{x, x')M-\x', 0) = M(x, x')M-\x', y)M~\y, 0) = M'^x, 0) , (51) 

one replaces the r.h.s. of eq. (^) by M^^{x, 0) and solves this equation with the standard 
method. Thus, the computational effort to determine M~^{x, 0) is just twice the effort for 
the quark propagator. 

We didn't comment yet on the structure of Cnqn in terms of quark propagators. As 
the difference to GNqN is only that the sum over yo is missing, we could write down the 
connected and disconnected parts of GjygN exactly in the form of eq's. (|49|), (^OD, and (|47|), 
albeit omitting the sum over yQ. We could also use the two step procedure outlined above 
for the connected part. However, in doing so one would have to apply the second step of 
this procedure for each time slice separately, increasing the numerical effort linearly with 
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the number of time slices. In order to avoid this increase one can use alternatively the 
'insertion' procedure invented by Martinelli and Sachrajda |6^. This method allows for an 
arbitrary variation of yo without additional inversions. 

We have seen already, that M~^ can be easily obtained from the solution of eq. (^Bf ) if 
the the source, i.e. it's r.h.s., is suitably modified. This concept can be extended to even 
more involved effective propagators. Martinelli and Sachrajda use the main part of the 
proton propagator as a source. Schematically, one has 



y 
M{ P(0)^:——^^PM } = 



P (0)'<:z Z^ P ^^^ ■ 

(52) 



The r.h.s. of this equation is given by a combination of quark propagators. Note that the 
timeslice Xq of the nucleon sink is fixed. CjyqN is obtained from the product of the solution 
to eq. (0), C{y, 0), and M~^{y; 0), for arbitrary yo. The expected ground state behavior, 
eq. (|39|) , can be identified as follows: One first chooses a large value of xq. Then one varies 
yo between and xq to obtain the plateau region. 

The disconnected parts of GiyqN and C^qjy have quite simple structures. All we need 
to know is the nucleon propagator and the 'quark loop', J2yTrM~^{y,y). However, if one 
would use eq. ( ^6l) to determine the trace, the computational effort would be A^^ x Nt times 
the cost for the calculation of a single quark propagator, a prohibitively large amount of 
computer time. 

The key to reduce this amount drastically can again be found in a suitable modification 
of the source term of eq. (|46|) . Suppose we would insert the source /, with ly = 1, to the 
r.h.s. . Then, the solution M~^ to this equation is given by 

M"i = Y,M-\x,y)Iy = Y.M-\x,y) . (53) 

y y 

Summation over x yields 

J2J2M-\x,y) = Y.TrM-\y,y)+ ^ M-\x,y) . (54) 

X y y x,y;x^y 

As the second term on the r.h.s. is not gauge invariant, one expects it to cancel out in the 
process of averaging over all gauge configurations. Note however, that this cancellation 
is mingled with the 'normal' statistical fluctuations of J2yTrM~^{y]y). Since one is in 
practice always limited to a finite sample of gauge configurations, a systematic bias due to 
a non-cancellation of the second term could be easily overlooked. Furthermore, one would 
expect an increase of the statistical uncertainty, since the noise of this term adds to the 
noise of J2yTrM~^{y,y). 



Thus, it is highly desirable to modify this 'volume source' method [^ such that it 
allows to remove J2x,y;x^yM~^{Xjy) on each single gauge configuration. This is achieved 
by the 'stochastic estimator' methods [^,^, which we will now explain briefly. 
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The major difference between the 'volume source' method and the stochastic estimator 
technique is that the volume source vector / is replaced by a random vector 77, which has 
the properties 

(?7i) = and {r]ir]j) = 6ij . (55) 

Note that the brackets denote the average over the distribution of random numbers on 
each gauge configuration. For each choice s of a random source rj^ one gets 



M-V = ^M-i(x,i/h 



(56) 



as a solution to the modified eq. (^Bf). The average over A^^; random sources then yields 



Ne 



Nf 



Y.Y.€M-\x,yX 



=1 x,y 



1 



Ne 

E 



Y.rilM-\y,yX+ ^ r^lM-\x,yX 



x,y;x^y 



■ (57) 



For large Ne the second term on the r.h.s. cancels according to eq. (|55|), and the first 
term reduces to the desired result TrM~^ on each single gauge configuration. The onset 
of asymptotics and the statistical accuracy of TrM~^ can be monitored within the course 
of sampling estimates. We emphasize that both these properties depend on the choice of 
the random number distribution. It has been shown in ref's. 0,|7^ that a (complex) Z2 

Unfortunately it turns out that 



distribution |]63[ is superior to a Gaussian distribution |K 
even then the number of estimates Ne has to be chosen of 0(300) for current lattice sizes. 
We mention finally that the stochastic estimator technique allows to compute single 
matrix elements of M~^. For example, using Z2 noise|^, the combination 



1 



Ne 






%M~ 



{xo,y)n 



1 



Ne 



— T 



yields M ^{xq, zq) in the limit A''^; -^ 00. 



M-\xo,zo)+J2v:oM-\x,yX 



(5J 



2.3 Lattice Results 

2.3.1 Quenched QCD 

The numerical determination of disconnected contributions by application of the methods 
outlined above requires both, high statistics of gauge configurations and, at least if one 
uses stochastic estimator techniques, a multiple of compute power (0(300)) compared to 
standard propagator calculations. Early, exploratory lattice simulations 0,^,|65| were 
therefore restricted to the calculation of the connected part of ct^at in the quenched ap- 
proximation. As a result, an astonishingly small value, a^°^ ~ 15MeV, was found in these 
calculations, suggesting that the disconnected contribution to the full pion-nucleon-sigma 
term should be approximately twice this number. 



^For Z2 noise, the property of the average {rjirji) = 1 holds for each estimate, i.e. rjiiji = 1. 
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Figure 3: Ref. ffB^: The ratio R^^^'(t) for the connected and disconnected contributions 
(u+d insertion) to the scalar density amplitude of a proton {P\uu + dd\P) , at a quark mass 
corresponding to mT,/mp = 0.604. 



Due to the advances in computer speed and algorithms, it is nowadays possible to 
investigate connected and disconnected parts in some detail. A high statistics simula- 
tion (300-400 gauge configurations, depending on the quark mass) has been performed by 
Fukugita et al. |^ at a lattice cutoff of a~^ ~ 1.5GeV. They used a slightly modified 
wall source technique for the disconnected contributions and analyzed all signals with the 
summation method, eq. (^Tj). 

To illustrate the statistical quality of their data we show in fig. |^ the ratio R^'^'^^{t) 
both for the connected and the disconnected contributions to cTttat. Note that, although 
the errors to the disconnected data are large, the expected linear rise of the signal can 
nevertheless be identified in both cases. 

A similar analysis, although with a significantly smaller statistical sample (24 gauge 
configurations for the connected and 50 for the disconnected part) has been done by Dong 
et al. 1^ at a larger cutoff, a^^ ~ 2GeV. They have applied the stochastic estimator 
technique, eq. (|F7|), with (complex) Z2 noise and 300 estimates per configuration. 

In fig.^ we show the results of both simulations as a function of (m.,^/mpY . 

The data for the connected insertion is consistent within statistical errors, whereas the 
estimates of ref. [^ for the disconnected contribution appear to be significantly lower 
than those of ||BB|. This might indicate a finite cutoff effect, but in view of the very limited 
sample of gauge configurations used in ref. [0 it is likely that the statistical errors have 
been underestimated substantially. 
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Figure 4: Connected and disconnected contributions to the scalar density amplitude of 
a proton, {P\uu + dd\P). The data has been taken from Fukugita et al. / |g^/ and from 
Dong et al. pl^. To compare the results of different simulations we have renormalized the 
amplitudes consistently with the tadpole improved renormalization factor Zs Wl 



Obviously, the mass dependence of the connected and the disconnected part is quite 
different. The disconnected part grows significantly with decreasing quark mass whereas 
the connected part is basically unaffected by the variation of mass. 

An extrapolation (see below) to the hght quark mass, ('^7r/'"^p)ph,j;s = 0.033, consoli- 
dates the presumption from the results of the early quenched calculations. One finds 



Tjquen 



{P\uu + dd\P) 



disc 



{P\uu + dd\P), 



2.23(52) 
1.8(1) 



ref. 
ref. 



[66] 




(59) 



and 



quenched 



(60) 



■ 43(7)MeV ref. [66] 
^ 50(3)MeV ref. H ' 

Within errors the quenched findings of a-j^N are consistent with the estimates extracted 
from experiment. The disconnected contribution appears to be much larger than the con- 
nected one. Thus, according to eq. (p^, one needs a strong breaking of flavor symmetry 
for the disconnected part to achieve a low value of y. 



Note that a^.^ is a renormahzation group invariant quantity. Thus, the lattice result, 
multiplied by the lattice scale yields directly the physical result. 

The matrix element and the quark mass individually are connected to their respective 
physical values by renormalization factors Zg = 1/Zm, which depend on the bare coupling 
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go and the lattice quark mass 

{P\uu + dd\PY^y' = Zs{P\uu + dd\PY'''' , mfy' = Zramf' . (61) 

We will discuss the methods to determine these factors in the next chapter. Here, we 
mention however, that the lattice result for the light quark mass used to obtain a-^^N cor- 
responds to a physical mass m/(2GeV)^/5 ~ 5MeV. 

We did not comment yet on the ansatz used for the extrapolation of the matrix elements 
to the light quark mass. Since both simulations have only three data points at their 
disposal, and in view of the fact that the statistical errors to the disconnected parts are 
large, one would not expect to be able to resolve reliably contributions from higher orders 
in the quark mass. Correspondingly, ref. |BB[ uses a linear ansatz. 



Inspired by chiral perturbation theory, ref. |^ extrapolates with a non-linear, non- 
analytic ansatz of the form C + D^Jm^. This, of course, leads to an enhancement of 
the results at the light quark mass. However, even in high precision calculations of the 
hadron spectrum |6^, with 0(800) gauge configurations and quark masses corresponding 



to {mT^/iTLpY = 0.16, one does not find numerical evidence for the presence of such non- 
analytic contributions. 

We have seen above, that the ratio of disconnected to connected contributions to o",rAr 
is astonishingly large in the quenched approximation. Assuming an approximate flavor 
symmetry of the disconnected parts, one would conclude that the y parameter, defined 



in eq. (piq) , takes the value y = 2/3. This is much higher than the phenomenological 



expectation, y ~ 0.2 — 0.4, and it would indicate that the strange quark contribution to 



the nucleon mass, ms{N\ss\N) is as large as ~ 400 MeV [IS]. 

To determine y and to check the assumption of flavor symmetry, one has to calculate 
the (purely disconnected) strange density matrix element of the proton. In the quenched 
approximation, this is done in two steps. 

First, one adjusts the lattice quark mass such that the physical strange meson spectrum 
{K,K*,(^) is reproducecQ. This determines the strange quark mass. In the second step 
one extrapolates the matrix element as a function of the valence quark mass of the proton 
to the light quark mass, while keeping the loop quark fixed at the strange mass. This then 
yields the lattice value of (A^|ss|A^). 

Using this procedure, ref. [^ finds 

y = 0.66(15) , (62) 

which is almost exactly the value one expects from the assumption of flavor symmetry 
of the disconnected contributions. The result appears plausible, since the quenched QCD 
vacuum is not sensitive to quark mass and quark flavor. 



^^Tlie value for the strange quark mass depends (O(10%))on the physical quantity {rriK or mx" ,m^) 
used to define it |42|. This however does not influence the results discussed below. 
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Note that y depends on the ratio Zs{ms)/Zs{mi). Since the change of Zg with the 
lattice quark mass is expected to be small in the light and strange quark region, this ratio 
should be close to 1. Using the tadpole improved perturbative renormalization scheme p9[ , 
which has been applied in ref. ||66|, one finds Zs{ms)/Zs{mi) ~ 1.06. 



It has been suggested by Lagae and Liu [^ that the lattice quark mass dependent 
part of the renormalization factor would be quite different for connected and disconnected 
amplitudes. Use of their result would induce a strong symmetry breaking in the discon- 
nected sector. The authors of ref. ||6^ have applied this modified renormalization procedure 
to their raw data. On top of this they have used the non-analytic ansatz for the chiral 
extrapolation discussed above. Both manipulations lower the value of y, and they find 
y = 0.36(3). 

We do however not agree with the suggestion of Lagae and Liu, at least in the case 
of scalar quark loops, for the following reason: Since the combination mq{P\qq\P) is a 
renormalization group invariant quantity, both parts of the amplitude, connected and dis- 
connected, will be renormalized by the same factor Zs = l/Zm, where Zm renormalizes 
the quark mass. Choosing different factors, Z|.°" 7^ Zf'^'^ would violate the renormalization 
group invariance. 

In summary, quenched lattice simulations find ct^at to be consistent with the phe- 
nomenological expectation. It turns out that the ratio of disconnected to connected con- 
tributions is close to 2 and that the fiavor symmetry of the disconnected part in nearly 
unbroken. This leads to a large value of the y parameter. It is therefore of utmost interest 
to study these aspects in full QCD lattice calculations. 

2.3.2 Full QCD 

We have seen above, c.f. fig. ||, that quenched simulations with a statistics of 0(300) 
configurations find statistical errors on disconnected contributions of ~ 30%. This is bad 
news if one plans to determine such contributions in full QCD, for two reasons. First, 
the creation of gauge configurations in full QCD enhances the computational demands by 
a factor of 0(1000). Thus, one needs a 'high end' supercomputer to generate statistical 
samples with comparable size. Secondly, a comparison of quenched and full QCD spec- 



trum results at equal statistics revealed |]70[, that the statistical fiuctuations on full QCD 



configurations might be much larger than in quenched QCD. Early, exploratory full QCD 



simulations l7T|,^] have therefore not been able to extract possible unquenching effects on 
a^N and y. 

Recently, the SESAM collaboration has generated 4 samples of 200 statistically inde- 
pendent gauge configurations with nj = 2 mass degenerate Wilson quarks, on a n^ x nj = 
16^ X 32 lattice, and a lattice cutoff a~^ ~ 2.3GeV. The samples correspond to different 



sea quark masses, with m^/nip = 0.833(3), 0.809(15), 0.758(11) and 0.686(11) |70| . 



The analysis [0,|74[ of these configurations indeed showed that the standard procedures. 



i.e. summation method combined with the wall source or the stochastic estimator technique 
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- fail in producing a reasonable signal for the disconnected contribution to Ct^n- The 
situation can be significantly improved however, if one uses PAM, together with a stochastic 
estimator technique with (complex) Z2 noise. 

We illustrate this result in figs. ^ and |^ , where we compare the signals for the ratio 
jlPAM ^ created using the Z2 stochastic estimator method with those for R^^^\ generated 
with the wall source and, alternatively, with Z2 stochastic estimator technique. 

Obviously, the former combination (/?^^*^+ Z2 stochastic estimator) is superior with 
respect to the statistical quality of the signal. However, the errors are even then of 0(30%). 





Figure 5: Disconnected contribution to a-^N- Comparison of wall source + summation 
method with Z2 stochastic estimator + PAM at a sea quark mass corresponding to rriT^/mp = 
0.833. The solid lines indicate the region of ground state behavior of the signal. 



Using the plateau method for the connected contributions and PAM + stochastic Z2 
noise for the disconnected part, the authors of ref. |]73| have determined the scalar am- 
plitudes of the proton at the 4 SESAM quark masses. The results are shown in fig. ^, 
together with the linear extrapolations to the light quark mass. 

At first sight the full QCD results look quite similar to those of the quenched analysis, 
c.f. fig. ^. To see this in some more detail we compare in tab. |I| the Uf = 2 results at the 
light quark mass with the findings of of ref. |^ . 



It turns out that, although the amplitudes are consistent within errors, the relative 
weight of disconnected and connected amplitudes is quite different. One finds 



R 



full 
dc 



1.26(57) 



(63) 



compared to -R^"*^" = 2.23(52). This might give a first indication for the influence of sea 
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Figure 6: Disconnected contribution to a^N- Summation method with PAM at a sea quark 
mass corresponding to m.,^/mp = 0.686. In both methods, the stochastic estimator technique 
with complex Z2 noise has been used. The solid lines indicate the region of ground state 
behavior of the signal. 





Zs{Cu + C,) 


Zs{'iD,) 


sum 


Zs^mia-\MeY] 


a"i[GeV] 


ref. HI 
quenched 


2.62(6) 


5.8(1.4) 


8.6(1.4) 


5.0(2) 


1.45(2) 


ref. jTSl 


3.11(25) 


3.9(1.8) 


7.0(2.0) 


2.7(2) 75| 


2.30(6) 



Table 1: Results at the light quark mass. For comparison all data has been renormalized 



by the tadpole improved perturbative renormalization factors Zs, Zm JH/ 



quarks on scalar density nucleon matrix elementsj^. As we mentioned above, the ratio Rdc 
is one of the two quantities which fix the value of the y-parameter. 



A consistency check to the sum of amphtudes can be done in full QCD by use of the 
Feynman-Hellmann theorem, eq. (|3^) . From a quadratic fit to the proton mass data, the 



""^^Some caution in the interpretation of this difference is in order since the results have been obtained 
at unequal lattice cutoffs. 
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Figure 7: Full QCD: Linear extrapolation of connected and disconnected amplitudes to 
the light quark mass, from ref. jfT^. Bursts indicate the results of the extrapolation. Cq 
denotes the connected , and Dq the disconnected contributions to (Jt^n- 



authors of ref. ||73] find 



Zs^{rni) = 9.2(3.9) . 



(64) 



This is in agreement with the result of the direct determination. Note, that the statistical 
error in eq. (0) is larger. 
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From tab. |1] one extracts a rather small value of the pion-nucleon-sigma term 

a^N = 18(5)MeV . (65) 

Clearly, this is almost exclusively due to the weight of the light quark, which is lower by 
nearly a factor 2 in the Uf = 2 case compared to the quenched result. We emphasize, that 
in both cases the standard lattice definition]^ of quark masses [^] has been used. 

It has been argued by R. Gupta |7^ that the perturbative renormalization would be 



inappropriate in full QCD, and that non-perturbative contributions due to the presence of 
sea quarks would enhance Z^ by approximately a factor of 2. This however does not affect 
the result for ct^at, as this is a renormalization group invariant quantity. 

Another explanation could come from the lattice definition of the quark mass. Due to 
finite cutoff effects, different definitions of physical quantities can lead to different answers 
on the lattice. Of course, the results have to agree in the continuum limit. 

It is known from quenched simulations, that the standard definition of the lattice quark 
mass, which is extracted from the discretized Hamilton density, and the PCAC definition, 
which uses the axial Ward identity fT^, do not agree at finite lattice cutoff. At a'^ ^ 2GeV, 
the (perturbatively renormalized) quark masses differ by about 30%. A scaling analysis 
revealed [^] that the cutoff dependence of the quark mass in the standard definition is 
much weaker than in the PCAC definition, thus favoring the former. 

The situation could be different in full QCD simulations. Indeed, the (preliminary) re- 
sults of the CP-PACS collaboration [^ indicate a stronger cutoff dependence of the quark 
mass in the standard definition. For a definite conclusion on that issue however more data 
at different cutoffs are needed. 

We now turn to the determination of the y parameter in full QCD. In principle, since 
one needs a strange quark loop contribution in this context, a fully consistent QCD lattice 
calculation would require both, light and strange quarks to be present in the quark sea. 
Such a n/ > 3 simulation with non-degenerate quarks will surely be one of the next large 
projects in lattice QCD, but is has not been performed yet. 

Within the framework of QCD simulations with degenerate light sea quarks, one has 
to rely on a semi-quenched analysis, where the strange quark has no counterpart in the sea 
of quarks. This at least accounts for the light quark mass dependence of the vacuum and 
is therefore more appropriate than a pure quenched analysis of y. 

The semi-quenched method is technically very similar to the quenched procedure out- 
lined above. The physical difference is that the mass of the light valence quark is always 
identified with the sea quark mass. 

Using this procedure, the SESAM collaboration finds |7^ 

Z5(P|2ss|P) =4.1(1.5) , (66) 



^^The standard definition of the bare quark mass is given in terms of hopping parameters k by iriqa = 
1/2(1/ K— 1/kc)- The critical parameter Kc is determined such that the pion mass vanishes for this choice 
of K. 
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which is in accordance with the result from the hght quark region, given in tab.^ This 
indicates, that the flavor symmetry of the disconnected contributions to the scalar proton 
amplitude is still intact, even in the presence of light sea quarks. 

Thus, only Rdc serves to reduce the value of the y parameter. Correspondingly, the 
result 

y"/=2 = 0.59(13) (67) 

differs only slightly from the quenched value. 

2.4 Discussion 

Let us shortly comprise the physics information that evolves from lattice QCD simulations 
discussed above. Taking the results at face value, the following picture arises: In quenched 
QCD, ajrN is consistent with the 'experimental' value. However, the major contribution 
(2/3) comes from the vacuum polarization (disconnected) part of the amplitude. As those 
parts turn out to be flavor invariant, the y parameter picks up a rather large value (2/3), 
in sharp contrast to the phenomenological expectation. 

In full QCD, with two degenerate sea quarks, the connected part is slightly enhanced 
whereas the disconnected amplitude is decreased. The sum of both contributions is still 
close to the quenched result. Because the light quark mass is found to be quite small, the 
full QCD lattice result differs from the 'experimental' value by about a factor of two. The 
flavor symmetry of the disconnected parts seems not to be hampered by the presence of 
light sea quarks. Thus, the y parameter decreases only slightly. 

Let us now discuss what could be done to consolidate this picture. First of all, the 
statistical errors to the disconnected contributions are too large. Although the stochastic 
estimator technique with complex Z2 noise in combination with PAM reduces these errors 
considerably, one is still confronted with uncertainties of 0(30%), on gauge conflguration 
samples of 0(200). 

In the case of quenched QCD it should be no major problem to reduce these errors by 
a factor of 10, simply by exploiting the computational power of currently available parallel 
supercomputers. 

For full QCD calculations the situation looks less promising. One can expect however, 
that the signal to noise ratio will improve on larger lattices, due to self averaging effects. 

Once one has succeeded in reducing the statistical errors to a few percent, a scaling 
analysis of cryrAf and y can be done. This will remove possible cutoff effects and lead 
to reliable continuum values for a-^^N- Since the y parameter depends on (a ratio of) 
renormalization factors it will be advantageous to use non-perturbative renormalization 
techniques ||78|,|7^ instead of the standard (tadpole improved) perturbative methods. This 



might weaken the cutoff dependence of y considerably. 

The last step in a full QCD determination of y will be to include sea quarks of dif- 
ferent masses into the generation of vacuum configurations. A lattice simulation with 
two degenerate light and two degenerate strange quarks should be feasible even with cur- 
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rent algorithmic and computational facilities, as it avoids the sign problem of the fermion 
determinant 
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3 The Flavor Singlet Axial Coupling of the Proton 



About a decade ago the European Muon Collaboration (EMC) ||T9l announced measure- 
ments of the first moment of the polarized proton structure function Qi, which implied an 
unexpectedly low value for the flavor singlet axial coupling of the proton G\ ~ 0, defined 
by 

s,G\ = {p,s\A'^{q)\p,s),2^o- (68) 

s^ is the proton spin vector, and A^ is the flavor singlet axial current 

K = ^A E ^7/.75g. (69) 

q=u,d,s 

Z\ denotes the renormalization constant of the singlet current in the continuum. 

In the simple parton model, G\ can be interpreted as the fraction of the proton spin 
carried by the quarks. Thus, on first sight, the EMC result indicated that the quarks would 
not contribute substantially to the proton spin: The confusion which followed from this 
naive interpretation was called the 'proton spin crisis'. 

On top of this, the EMC result stands in contradiction to the Ellis- Jaffe sum rule [^ 



which can be derived in the OZI (Zweig rule) limit of QCD. This indicated already that the 
proton spin problem might be deeply connected to the quantum field theoretical properties 
of QCD. 



Since then, in a tremendous effort, both on the experimental pO|-M and on the theo 



retical side |81], the EMC result has been refined and a clear qualitative understanding of 



the phenomenon could be achieved. 

The latest analysis of experimental results [^, including all proton, deuteron and 
neutron data, yields]^ 

G\ = 0.29 ±0.06, (70) 

which is significantly below the Ellis- Jaffe prediction of {Gr'j^^'^ = 0.58. 

From the theoretical investigations it appears that it is the influence of the U{1) axial 
anomaly, which leads to the small value of G\. This immediately removes the 'proton spin 
crisis', as an interpretation in terms of the simple parton model is no longer valid. We 
mentioned in the introduction that the axial anomaly is closely connected to the topolog- 
ical properties of the QCD vacuum. Thus, the anomalous value of G\ is one of the rare 
opportunities for a quantitative study of topological aspects. 

The lattice approach provides (at least) two different methods for the calculation of 



G\. 



The first one, which we will call the direct method, is to evaluate the r.h.s. of eq. 
1). Since this is a flavor singlet matrix element, one expects contributions from both 



-"This is at a renormalization scale of Qq — 5GeV^. 



S. Giisken, Flavor singlet phenomena in lattice QCD 



31 




N(p) 



X(p') 



Figure 8: Inclusive polarized DIS scattering. 



the connected and the disconnected parts, depicted in fig. |I[ If it turns out, that the 
smallness of G\ originates from sizeable (negative) disconnected contributions, one could 
immediately clarify why and to what extent the Ellis- Jaffe sum rule fails: Disconnected 
diagrams are excluded in the OZI limit of QCD. 

The second method, which we will denote as the topological method, focuses on the 
topological aspects - charge and susceptibility - of QCD. Here, one calculates the cor- 
relation of the proton propagator and the topological charge density. A success of this 
approach with respect to a small value of G\ would encourage for a detailed, quantitative 
investigation of the connection between the axial anomaly, QCD topology and instantons. 

Both methods will be explained below in detail. 



3.1 Preparation 

Let us briefiy recall how G\ is connected to the polarized deep inelastic scattering (DIS) 
of leptons and protons, shown schematically in fig.^]. 

The differential cross section for lepton proton DIS in the laboratory frame is given by 



dVLdE' 






(71) 



where M is the proton mass, E and E' are the energies of the incoming and outgoing 
lepton, and L^^, and W^^i, denote the leptonic and hadronic tensors, respectively. W^i, can 
be parameterized by two spin independent structure functions, Fi and F2, and by two spin 
sensitive functions, g^ and (/f 
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s denotes the polarization vector of the proton. 

The polarized structure functions gf and f^f can be extracted experimentally from the 
measurement of the asymmetries 



All 






Ai 



da^^ — da'^ 
da^^ + da^~ 



(73) 



where the arrows indicate the polarizations of the lepton (first arrow) and the proton 
(second arrow): up and down arrows stand for longitudinal polarizations, and the right 
arrow denotes a proton spin direction transverse to the lepton momentum and towards the 
direction of the scattered lepton. 

On the theoretical side, gi is given by the proton matrix element of two electromagnetic 
currents carrying large space like momentum. Its first moment 



niQ')^ f'dxg',{x,Q'), Q' = -q\ x=^ 
Jo 2va 



(74) 



can be derived using the operator product expansion (OPE) [^. Including terms up to 
twist two one obtains 



JP{q)r{-q)^"^^2e'^--'^^ 


cricks) [ai + i=4) + IcUaMl 


and correspondingly 


r?(Q^) = ^c( 


'' {g\ + Ig\) + lc^G\ . 



(75) 



(76) 



The Wilson coefficients C^^ and Cf are both known to 0{a^) |^. The non-singlet 
form factors G\, G\ and the singlet coupling G\ are given by the forward proton matrix 
elements of the (renormalized) axial currents 



{p,s\A^\p,s) 

{p,s\Al\p,s) 
{p,s\Al\p,s) 






(77) 

(78) 
(79) 
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with 

4 = 2r^(^7^75W + rf"7/.75C?-2s7^75s) (80) 

< = ^i E ^7M75g- 

q=u,d,s 

Note that we have introduced different renormahzation factors for the non-singlet (Z^^) 
and the singlet (Z^) axial currents. This is necessary, since the conservation of A^'^ is only 
violated by the quark mass, whereas the singlet current is not conserved in the limit of 
vanishing quark mass, due to the presence of the axial anomaly. 

The non-singlet couplings can be extracted from neutron and hyperon (3 decays. One 
findsQ BH 

r^ — ±.Z,U ( U _1_ u.uuou , u^ 



Gl = 1.2670 ± 0.0035 , Gl = 0.579 ± 0.025 . (811 



Thus, the flavor singlet coupling G\ is determined by the measurement of Ff . 

In the OZI limit of QCD, where (strange) quark loops do not contribute, eq's. ( ffSD 
and ([79|) imply G\ = G\, in sharp contrast to the experimental result and to the fact 
that the Q^ dependence of G\ is different from that of G\. The QCD parton model can 
qualitatively repair this failure by the assumption of a large gluonic contribution to Qi. It 



turns out that singlet quark and gluon distributions mix under renormahzation |8^. In 
a particular renormahzation scheme, the Adler-Bardeen scheme, it is possible to shift the 
anomalous Q^ dependence completely to the gluonic part, leaving intact the naive inter- 
pretation of the singlet quark distributions as the fraction of the proton spin carried by the 
quarks. However, the QCD parton model is unable to make any quantitative prediction 
for G\. 

In this situation, a non-perturbative method, which has direct access to the singlet 
proton amplitude, c.f. eq. (|79D , is badly needed. A lattice computation of this amplitude 
can reveal whether, and to what extend, QZI violating contributions determine the value of 
G\. In the language of fig.|l], the latter correspond to the disconnected parts of the singlet 
proton amplitude, whereas the connected contributions should still obey the Ellis- Jaffe 
sum rule. 

Before we explain how such a direct calculation can be performed on the lattice we will 
discuss shortly the second approach to G\ which exploits the close connection between 
flavor singlet axial current, axial anomaly and topological charge. 

The axial vector Ward identity 

^M< = 2^ E ^<l.Q^l5q^ + 2nf^trG^,G^'' , (82) 

i=l,nf 



^^G\ has been determined from hyperon decays under the assumption of SU{'S) flavor symmetry. Thus, 
the uncertainty on this quantity might be underestimated. 
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with Hf being the number of quark flavors and G{G) the (dual) tensor of gluon field 
strength, states that the axial vector flavor singlet current is not conserved in the chiral 
limit due to existence of an anomalous pure vacuum contribution trGfj^^G^'^, called the 
axial anomaly. 



The anomaly is deeply connected to the topological structure of the vacuum |^^. In 
particular, the topological charge Q (topological index, Pontryagin index), which labels 
the different topological sectors of the vacuum, can be expressed by 

Q = -j^ld'xtrG,,G^''' . (83) 

Thus, the axial anomaly is just the topological (charge) density Q{x) = trG^^G^^ of the 
vacuum. 

We can now reformulate eq. ([79| ) in terms of Q. In the chiral limit (m^ = 0), a Fourier 
transformation of eq. ( p^ ) yields 

iq^.{v\s\Al{q)\p,s) = 2nf^{p\s\Q{q)\p,s) , (84) 

where p,p' are the 4-momenta of the incoming and outgoing proton, and q = p — p^ . If 
we choose the proton polarization s to be (anti)parallel to the direction of the momentum 
transfer g, 

we obtain in the rest frame of the proton at zero momentum transfer 

G\ = s,{p,s\Al\p,s) = \im2nf-^^^{p',s\Q\p,s) . (86) 

g'-»o 4:71 q ■ s 

Thus, one arrives at a topological expression for the singlet axial coupling. We call the 
lattice procedure to determine the disconnected matrix element on the r.h.s. the topological 
method. Details will be given below. Let us mention however, that this method has a 



serious shortcoming in the quenched approximation [Q. Because of the fact that the rj' 
meson becomes massless in the chiral limit of quenched QCD, the desired matrix element 
{p',s\Q\p,s) diverges in this limit. Therefore, such calculations can be done only in full 
QCD. 

An alternative, less fundamental, expression for G\ in terms of topological quantities 
can be found using the composite operator propagator- vertex (CVP) method |jT5|, |16|, ^ . 



Here one uses the OPE to decompose the singlet part of F^, eq. (0), into composite 
operator propagators {0\TOiOj\0) at zero momentum and vertex functions F 

in)s^n, = l^^^fCtic^s) [(0|rQQ|0)FQ,, + (0|rQ$5|0)F^,,,] , (87) 

with $5 = J2i=i,nf qil^qi- The topological susceptibility 

X{k^) = Jdxe'^''t{0\TQ{x)Q{0)\0) (88) 
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vanishes in the chiral hmit of full QCD at zero momentum transfer. The remaining prop- 
agator on the r.h.s. of eq. ( ^7]) is given by the square root of the first moment x' of the 
topological susceptibility [1^. Thus, if the vertex function F^^pp and the first moment of 
X would be known, one could calculate the singlet axial vector coupling. Unfortunately, 
the vertex function cannot be calculated directly from first principles, making the CPV 
approach somewhat less rigorous than the topological method. However, it can be estimated 
e.g. in the QZI limit of QCD. The determination of x' is a typical task for non-perturbative 



methods like lattice QCD or QCD spectral sum rules (QSSOR). The authors of ref. [16| 



have used the latter. They founcQ y X' i^) / Q2=ioGeV^ = 23.2 ± 2.4 MeV and, correspond- 
ingly G\ = 0.35 ±0.05, consistent with the result of the latest experimental determination, 
quoted in eq. (|70D . It will be interesting to see, whether this result can be consolidated by 
a full QCD lattice computation. 

Finally, let us comment on the apparent incompatibility in the physical interpretation 
of the two approaches. From a success of the direct method one would conclude that 
the correlation of vacuum quark loops with the proton is responsible for the small value 
of G\, whereas one expects in the topological method that topological excitations of the 
vacuum decrease the flavor singlet axial coupling. We want to emphasize that these are 
not different mechanisms: Topologically non trivial vacuum configurations influence of 
course the correlation of quark loops with the proton, since it is exclusively mediated by 
the vacuum. From that point of view, one would consider the (disconnected) quark loop 
contributions to G\ as being induced by the topological properties of the vacuum. 



3.2 Direct Method 

According to eq. ( ffOD we need to compute the matrix element {p, s\q'y^'y^\p, s). The tech- 
niques which are employed for this task are quite similar to those introduced in section [2.1| . 
Again, one exploits ratios of 3- and 2-point correlation functions at large time distances, 
namely 



XQ- 



' A+(Ar|g7^75g|iV)xo, (89) 



yPLAf^ „. ,,^ _ ^NqNi^O^Vo^fJ') 



Ra ixo,yo,fj.) = — ^ ^ (IZ(g7M75g)(y,l/o)) 

'~JNN[Xo) g 



Xq^OO 



{N\q-i^^^q\N) , < yo < a;o , (90) 



^^A very preliminary result of an exploratory full QCD lattice simulation with Kogut-Susskind fcrmions 
has been reported by the authors of ref. Q: ^x'(O) — 16(3)MeV . 
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and 

XQ-Ax'g 

i?^^^(xo,Axo,Ax[„/i)= Y. i?A^^(a;o,l/o,/i). (91) 

yo=Axo 

The axial 3-point functions Gj^^^ and Cf^^^ are defined as 

Ci^,7v(a;o,Z/o,/i) = EE(^^(0.0)(g7M75g)(y,l/o)iV(x,xo)) , (92) 

X y 

GiqNixo,f^) = E<^ivq^(^o,2/o,/i) • (93) 



yo 



One might suspect that the apphcation of the summation m,ethod is not justified, since 
a term ~ ^7^755' is not present in the QCD action Sqcd, and thus the Feynman-Hellmann 
theorem cannot be used directly to derive eq. (|89|) . However, if one starts from a slightly 
modified action 

SqcdW = Sqcd + Xj2iQ'^Q)iy) , T = 7^,75,7^75, ... (94) 

y 

and repeats the steps which lead to eq. ( pi]) one arrives at 

^A (a;o, A) = Q ,^ y. - {2^{Q^Q){y)) = ^(A) + — Q^ — xo . (95) 

The Feynman-Hellmann theorem, see eq. (0), yields 

^^^ = (iV(A)|gTg|iV(A)). (96) 

Finally, in the limit A ^ 0, one obtains eq. (|89|). 



The numerical methods used to calculate the connected contributions to R^^ and 



R^^^ are identical with those described in section Bl2 



To determine the disconnected contributions, one has to evaluate 

iC^.Nf'" = {Fg^. (0, x) E Tri^,l,M-\y, y,- y, yo))) (97) 

y 

and 

{G^^Y-'^ = Y.iC^.NY'"' ■ (98) 

yo 



In principle, the standard estimator techniques, described in section |2]^, can be applied. 
However, these methods become inefficient if one computes off-diagonal elements of M~^ 
in spin space, which occur for example in the combinations Tr(7i 275M~^). In these cases, 
large diagonal contributions add to the 'noise term' on the r.h.s. of eq. (|58|), and one would 
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have to enhance the number of estimates or gauge configurations drastically to suppress 
these terms. 



Two alternative ways to cope with this problem have been proposed in ref s. |73, 91 



and ref. [^. The first method simply circumvents the problem by splitting the stochastic 
inversion, eq. (|56|), into four pieces, each of them corresponding to one spin component. It 
turns out that, with an equal computational effort, this spin explicit method reduces the 
statistical error on the estimate of Tr(7i^275Af "^) by about 40 - 50%. 

The second technique reduces the fluctuations from large diagonal contributions by 
subtracting the estimate for the on-diagonals 

1 Ne 

E[M-\xo,xo)] = ^T.T.v:,M-\xo,y)Vy (99) 



Ne 



s=l y 



from the estimate E[M "^{xq^zq)] for the off-diagonals, given by eq. ([581). This yields an 
improved estimate 



1 Ne 

E[M-\x,,z,)r^^ = E[M-\x„Zo)] - E[M-\xo,Xo)]—Y.V:Xo ■ (100) 



We emphasize that this diagonal improved method works equally well as the spin explicit 
method. 

3.2.1 Renormalization 

The axial current A^ is not conserved in the continuum due to the presence of flnite quark 
masses and to the axial anomaly. Thus it gets renormalized by quantum fluctuations. The 
effects of the two chiral symmetry breaking mechanisms on the renormalization of A^ are 
however quite different. 

The quark masses induce only a soft symmetry breaking, i.e. the dimension of the 
relevant operator is less than 4, and the renormalization factor Z^^ is cutoff independent 
and flnite. In particular, one obtains Z^^ = 1 in the limit of vanishing quark mass. Non- 
singlet axial currents are renormalized by Z^^ alone, as the axial anomaly is absent in this 
case. 

The symmetry breaking caused by the axial anomaly is not soft (dimension=4), and 
An picks up an additional renormalization which diverges logarithmically if the cutoff is 
removed |9^. This leads to a non-zero anomalous dimension which induces a renormaliza- 



tion point dependence on Z^ |^4[. Fortunately, the renormalization point dependence of 



G\ is weak. Increasing the renormalization scale /i^ from lOGeV^ to inflnity decreases the 
value of G\ by about 10% H. 

The situation is more complicated in lattice QCD. Lattice discretizations of the fermion 
action, which avoid the fermion doubling problem, introduce additional chiral symmetry 
breaking terms. It has been shown by Karsten and Smit |^ for the Wilson type of 



discretization that these terms are essential in reproducing the correct form of the axial 
anomaly in the continuum limit. However, the renormalization of the (axial) currents 
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becomes non-trivial. Although the lattice value, Z^^ , is still finite, the condition Z^^ = 1 
at TUg = is no longer valid. 

Much progress has been achieved in the determination of the lattice values of Z^^ , 
which depend on the lattice coupling and quark mass. We will comment on the different 
methods below. 

Not much is known about the lattice renormalization factor Z^ of the singlet axial cur- 
rent. In first order lattice perturbation theory, the logarithmically divergent contribution 
from axial anomaly is absent, and Z^ = Z^^ |31|. Unfortunately, a two loop or even a 
non-perturbative determination of Z^^ has not been performed yet. In view of the large 
statistical uncertainties in the disconnected parts of fiavor singlet axial current matrix el- 
ements (see below) it appears to be justified to work with the assumption Z^ = Z^^ . 
However, one should keep in mind that a possible improvement of the statistical quality 
would have to go hand in hand with a more precise determination of Z^. 

The early estimates of Z^^ were based on first order lattice perturbation theory calcu- 
lations, where the bare lattice coupling was taken as the expansion parameter |31|,^. It 
was demonstrated however through several examples |5^ that lattice perturbation theory 



in the bare coupling qq is ill behaved, and that large higher order corrections can introduce 
considerable uncertainties on the renormalized matrix elements. 

The convergence of the perturbative expansion can be improved, if one uses a renor- 
malized coupling rather than qq , defined by a proper renormalization condition p9|,p6 
Defining an improved coupling ay by the static potential V 



Viq) = -^-P^ , Cf = -, (101) 



CfAnaviq) 4 

Lepage and Mackenzie [^ have worked out an improved perturbative expansion on the 
lattice, which in many cases exhibits a faster convergence than the expansion in qq. More- 
over they demonstrated that the improvement is mainly due to the inclusion of tadpole-like 
diagrams, which, if ignored, induce a slow convergence of the expansion. Most of the recent 
lattice simulations have used this 'tadpole improved' perturbation theory in first order to 
renormalize their results. We emphasize however that even with this method higher order 
contributions can be sizable. Thus, a scaling analysis of the renormalized results is still 
mandatory. 

To overcome this problem, non-perturbative renormalization techniques have been de- 
veloped, which avoid the use of lattice perturbation theory ||97 Hl00|| • For the determination 
of Z^^ two different procedures, the Ward identity method (WI) |^ and the regularization 
scheme independent method (RI) [^ are equally appropriate. The WI method exploits 
the fact that the Ward identities on the lattice have the same structure as in the contin- 
uum. This defines a set of non-linear equations for the renormalization factors, which can 
be solved numerically on a given lattice. 

Clearly, the WI method is restricted to those operators where Ward identities do exist 
for. The RI method is more general in that respect and can be applied in principle to 
arbitrary composite operators. In this scheme the renormalization condition for an operator 
O is imposed on the corresponding amputated Green function, computed between off shell 
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quark states of large virtuality. This defines a renormalized operator O^, and thus a 
renormahzation constant Zq, which is independent of the regularization scheme ||101|| . A 
disadvantage of the RI method is that one has to work in a fixed gauge. This makes it 



necessary to study the influence of lattice Gribov copies on the results ||102 |. 

Both non-perturbative methods imply a numerical lattice determination of several cor- 
relation functions. Thus, uncertainties due to finite size effects, statistical accuracy and, 
in the case of the WI method, ground state identification may still be present. 



One may ask of course to what extent the value of Z^^ depends on the method used 
for its determination. To illustrate this point we show0 in tab.|] the results for Z^^ from 
the different methods, calculated with the Wilson action in the quenched approximation 
and at several lattice couplings /3 = G/g^. 



p 


a-i GeV 103| 


SPT 


TIPT 


WI 98 


RI 103| 


6.0 
6.2 

6.4 


2.258(50) 

2.993(94) 

4.149(161) 


0.87 
0.87 
0.87 


0.78 
0.79 
0.80 


0.85(7) 


0.81(1) 
0.81(1) 
0.82(2) 



Table 2: Results for Z^ in the chiral limit at a renormalization scale fi = a~ in the 
standard Wilson discretization. Denotation: SPT: first order standard (bare) perturbation 
theory; TIPT: first order tadpole improved perturbation theory; WI: non-perturbative Ward 
identity method; RI: non-perturbative regularization scheme independent method. 



Obviously, the effect of variation of methods is not dramatic (5 — 10%), compared 
to the large statistical uncertainties on the disconnected part of the axial flavor singlet 
current. Given that, together with the fact that remormalization effects on G\ are small 
in the continuum, one can hope, for lack of better wisdom, that the one loop result for the 
renormalization of this current, Z^ = Z^^ , is a currently tolerable approximation. 



3.3 Topological Method 

To determine G\ from eq. (|8^) by a lattice calculation one needs to define a discretized 
lattice operator for the topological charge density Q^ which agrees with the continuum 
definition of Q in the limit a — ;> 0, e.g. 



a-*0 



Q^(x)^^a^Q(x)+0(a^ 



(102) 



23 



This table has been originally compiled by the authors of rcf . [|103[ . 
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A straightforward way to define Q^ starts from the discretized version of the field strength 
tensor Gj^^, 



- p,,^(x) - P-u,-^ix) - P-^^A^) - P^MA] , (103) 

where P^^,^(a;) is the product of gauge parallel transporters U{x) around an elementary 
plaquette 

Pfi,u{x) = U^{x)Uy{x + afl)U^^{x + afi + aP)U^y{x + az>) . (104) 

From that it follows for the topological charge density ||104|| 



2^(^) = i E W<x^^ (P,,.(x)P,.(x)) , (105) 

^ fiupa 

where if^^ptj is the Levi-Civita tensor for /i, u, p,a > 0, and e^^pa^ = —e^^^p^, etc. . 

To relate the lattice definition of the topological charge density to the respective contin- 
uum quantity one needs to renormalize Q^. In principle, this is a highly non-trivial task, 
since Q mixes with dpAp and with the pseudo scalar current J2ili ^T^iQil^Qi under renor- 



malization |T^. It has been shown however in first order lattice perturbation theory | |105 | 



that the mixing is quite small (< 5%) for bare lattice couplings go < 1. Thus one can 
concentrate on the determination of one multiplicative renormalization constant Zq, 

Q^ZqQ'^. (106) 



Unfortunately, it turns out that the one loop correction to Zq is large ||105| , |106|] , making a 



non-perturbative determination highly desirable. Such a non-perturbative renormalization 
method has been developed by the authors of ref's. ||107| , |108|| . Roughly speaking it consists 



of preparing a classical gauge configuration of given topological charge and then to add 
gauge fluctuations by a Monte Carlo process. A comparison of the integer valued 'classical' 
topological charge with the one extracted from the 'fluctuating' gauge configurations allows 
for the determination of Zq. 

On the other hand one can try to improve the lattice operator Q^, eq. ( |105|) , such 



that the higher order corrections to Zq are made small or might even vanish. The link 
smearing technique | 109 |, which supplies a non-local version of eq. (|105|) , reduces the higher 



order corrections substantially. However, the 'ultima ratio' with respect to renormalization 
properties is the so called geometric definition of Q^, given by M. Liischer ||110|| . Here, 



Q'-' is constructed such that the corresponding topological charge takes only integer values 
even in the case of 'fluctuating' gauge configurations. It follows that Zq = 1, i.e. higher 
order corrections are absent. 

A comparison of the different definitions of topological charge on the lattice has been 



given by the authors of ref. [Ill] in the quenched approximation. It turns out that, us- 



ing the non-perturbative renormalization procedure described above, all definitions lead 
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to consistent results. However, the geometrical definition and the definition of eq. ( |105| ) 
suffer from large statistical noise. Thus, for practical simulations, one should use^ the link 
smearing technique ||109|| to calculate Q^. 



Having succeeded in defining suitable lattice operators Q^ there remains still a poten- 
tially severe problem for full QCD lattice applications^ to extract reliable results from 
numerical simulations. It has been obtained in ref. ||1 14|| for the staggered formulation of 
lattice fermions that the standard updating algorithm, the Hybrid Monte Carlo, tends to 
loose its ability to move between different sectors of topological charge if the lattice quark 
mass is chosen at or below ma ~ 0.01, which corresponds to rriT^/mp ~ 0.57. Thus, in this 
region of quark mass, the results might be biased by a non-Gaussian distribution of gauge 
configurations with respect to topology. The authors of ref. [[115|] have performed a similar 
investigation using dynamical Wilson fermions with quark masses down to a value corre- 
sponding to rriT^/mp ~ 0.69, finding no significant decrease of the Hybrid Monte Carlo's 



mobility. At niT^/mp ~ 0.56 ref. [ 115 | obtains still transitions between different topological 
sectors, but with a reduced frequency. This indicates that either the number of Hybrid 
Monte Carlo updates has to be increased at small quark massQ or one has to find an 
alternative updating algorithm which does not suffer from this kind of immobility. 

Assuming that we have successfully created a statistically representative sample of 
vacuum configurations we are now prepared to evaluate the r.h.s. of eq. (|8^) on the lattice. 
Obviously, two steps are necessary for this task: First one has to determine the matrix 
element (p', s\Q'-'{q)\p, s) for several momenta p,p' and secondly one has to extrapolate the 
results to the point q = p — p' = 0. Concerning the second step we emphasize that it is 
advantageous to work in the brick-wall system 

\p\ = \p'\ , go = 0. (108) 



Within this setting, the (unwanted) matrix element qQ{p',s\AQ\p,s) is absent in eq. ( plf ) 
even for qy^O. Thus, there is no need to remove it by the extrapolation g*— > in eq. (PB|). 
To extract the matrix element {p',s\Q'-'{q)\p,s) one starts from the 3-point correlation 
function 

A(y,r,x,t) = (pt(o,0)T,Q(y,r)P(f,t)) , (109) 



^ We did not discuss the 'fcrmionic' definition of the topological charge density 

Q^ = m, (Tr (75 A/- 1)), (107) 

where m„ is the quark mass and M~^ the quark propagator, which has been proposed by Smit and 



Vink |112|, inspired by the Atiyah-Singer index theorem. Although stochastic estimator techniques enable 
for an accurate calculation |112] of the r.h.s. of eq.(|l07|), the renormalization properties of have not been 



studied yet in detail. 

^^ Quenched QCD does not suffer from this problem. 

^®Ref. [115| estimates an increase by a factor of four when decreasing the quark mass from rriT^/mp ~ 0.69 



to m-^/mp ~ 0.56. 
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where P,P^ are interpolating proton operators and Tg is the spin projection operator. The 
expectation value A can be determined on the lattice by calculating the correlation of the 
(spin projected) proton propagator with the topological charge density. 

In Fourier space one obtains for given final momentum of the proton p and momentum 
transfei0 q 

Mq. r, p, t) = E A{y, r, f , t)e^^V^^^ . (110) 

At large time distances t — r and r the interpolating proton operator projects only on the 
proton ground state. This yields 

i(g,r,p,t)*-^°° (111) 



\P^\Pip',s)){Pip,s)\P\0){Pip',s)\Q\P{p,s)) X e 



-E{p){t^T)-E{p')T 



The proton amplitudes (0|P"'"|P(p', s)) and (P(p, s)|P|0) are determined from the asymp- 
totic time dependence of the proton propagator Ap(p, t), see eq. (pS)). Using the condition 
\p\ = \p'\ (see above) one finally gets 

{Pip',s)\Q\Pip,s)) = ^ lip... ^i^^p^'f • (112) 



-T,T — *00 



From a computational point of view one might think that the topological method is 
superior to the direct method since it circumvents the very costly and statistically noisy 
calculation of disconnected (fermionic) contributions with stochastic estimator methods. 
However it has some serious disadvantages which might be quite difficult to overcome. 
First of all, eq. (plD, which connects the axial current matrix element with the topological 
charge density, is valid only in the chiral limit. Apart from the fact that this makes an 
extrapolation in ruq necessaryp^ one has to neglect symmetry breaking effects according to 
'^M.d 7^ 0. Secondly one has to calculate A and Ap at several momenta p, q. The signals 
from momentum nonzero correlation functions are however very noisy on the lattice and one 
needs high statistics and well suited guesses for the wavefunctions to extract the required 
ground state information. On top of this the extrapolation in q might have to span a 
quite long range. The smallest nonzero momentum on a lattice with periodic boundary 
conditions is given by 

271" _i , , 

Pmin = — X a . (113) 

With typical lattice parameters A^^ ^ 20 and a~^ ^ 2GeV one gets the rather large value 
Pmin ^ eOOMeV. 

Last nor least, as we mentioned above, the application of the topological method is not 
justified in the quenched approximation [p3]. In full QCD one has to fight with a possible 



^ Note that the initial proton momentum ff is imphcitly fixed by p* — p — q. 

^^The fermion matrix approaches a zero mode at niq = 0. Thus, the quark propagator diverges at this 
point. A direct lattice simulation at ruq = is therefore excluded. 
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stiffness of tlie updating algorithm at small sea quark masses with respect to tunneling 
between different topological sectors]^. 

Exploratory lattice studies have been performed in quenched |^, |116|| and in full 
|117| , |118|| QCD to check on the feasibility of this method. It turned out that the statistical 



quality of the momentum nonzero correlation functions and, in full QCD, the stiffness of 
the updating algorithm with staggered fermions are the major hindrances in extracting 
reliable results. Thus, being optimistic, one can hope that the topological method might 
become useful with high statistics full QCD simulations using Wilson-like fermions. Such 
simulations are being performed currently by the CP-PACS |7^ and T^L ||119|| collabora- 
tions. 



3.4 Lattice Results 

3.4.1 Quenched QCD 

The first steps towards a lattice determination of G\ have been done several years ago 
by calculating the connected contributions to the matrix elements {P\uu\P) and {P\dd\P) 
|120| , |121|| . Although the statistical accuracy of these simulations has been rather low 



(less than 40 gauge configurations have been used), the results have been quite encouraging 
as they yielded estimates for the non-singlet axial couplings G\ and G\ in the ball-park 
of the experimental findings, eq. (pi|). 

Recently, a more precise estimate of G\ in the quenched approximation has been cal- 
culated by the QCDSF collaboration [|122|| . They used the 0{a) improved Sheikoleslami- 
Wohlert action with a non-perturbatively determined [^] 'clover' coefficient Cg^, and non- 
perturbative estimates for the axial vector renormalization constant Z^^ . The simulation 
has been performed at 3 different values of the bare coupling, f3 = 6.0, 6.2, 6.4, which corre- 
spond to lattice cutoffs in the range a~^ ~ 2 — 3.5Gev. The size of their statistical samples 
varies from 0(500) gauge configurations at /3 = 6.0 to O(IOO) at /3 = 6.4. Fig. | displays 
the results for G\ as a function of the (square of) the lattice spacing a in units of the 
string tension y/K. It turns out that the cutoff effects in this range of the bare coupling 
are small (0(10%)). A linear extrapolation^ to the continuum limit, a = 0, yields a value 
for 0^4 which is slightly lower, but, within errors, still consistent with the experimental 
result. Thus, one can expect from quenched lattice simulations to find a value of G\, which 
is close to experiment. Turning this conclusion around, the deviation of a lattice estimate 
of G\ at a finite cutoff a~^ from the experimental result can serve as an indicator for 
systematic effects in the lattice calculations of the flavor non-singlet couplings. For want 
of something more precise, this deviation may be used also as an (lower limit) estimate for 
cutoff effects present in lattice determinations of G\. 



^^Of course this problem affects also the direct method. 

•'"Using an 0{a) improved action one expects the cutoff effects to appear at 0{a^). Thus, a linear 
extrapolation of 0^(0^) is justified. 
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Figure 9: QCDSF collaboration \122^] : Continuum extrapolation of G\{a^). The star 
indicates the position of the experimental result. 



In quenched QCD, two such simulations have been performed up to now with the stan- 
dard Wilson action, one by Fukugita et al. ||123|| and one by Dong et al. ||124|| . The former 



has been done with accurate statistics (0(300) gauge configurations), but at a somewhat 
large coupling, (3 = 5.7, which corresponds to a~^ ~ l.SGeV. The latter has operated 
closer to the continuum, ai (3 = 6.0, but with a very low statistics (24 configurations), 
which might lead to a substantial under-estimation of the errors. 

Using G\ as an indicator for the size of cutoff effects, one finds in both cases that the 
uncertainty of their lattice estimates for G\ is dominated by the large statistical noise 



found in the disconnected contributions. Before we discuss the results of ref's. [ p.23| , |124 | 
in somewhat more detail we will illustrate briefly why one indeed expects considerable 
statistical fluctuations in those parts. 



We have seen in chapter ^]^ that the disconnected contribution to ct^at is much more 
noisy than the connected one. This is because both, the noise from the estimator method 
and from the inherent quantum fluctuations add up on top of the disconnected signal. From 
a numerical point of view the situation is even more disadvantageous for the disconnected 
contributions to G\. This is illustrated in fig. |10|. According to eq. (^) the disconnected 



part of cr^AT is given by the correlator of the nucleon propagator and the sum over quark 
loops with Dirac components (/i, /x), fi = 1,2,3,4. Since the contributions are of similar 
size for each of the Dirac components, the summation helps to stabilize the signal. The 
disconnected part of G\, see eq. (PB[), however is determined by the difference of two such 
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(1,1) (2,2) (3,3) (4,4) 

[0^0]> - <- <!> - [0^oj> 




(b) 



(2,2) (3,3) 



Figure 10: (a) Disconnected contribution to a-^N', (b) Disconnected contribution to 

{P\q^3^,q\P). 



noisy correlators. 

In fig. |ll] we display the signals for connected and disconnected contributions to G\ 
as found by the authors of ref. ||123|| , at a quark mass corresponding to m.,^/mp ~ 0.745. 
The connected parts exhibit a clear linear slope, see eq. (|^), whereas the disconnected 
contribution is very noisy and disappears at t = 12. 

Using the results of linear fits to the signals at several quark masses one can perform 
an extrapolation to the chiral limit. The quark mass dependence of connected and discon- 
nected contributions is shown in fig. |12|. Note that the amplitudes have been multiplied 
by the 1st order tadpole improved renormalization constant Z^ = Z^^ . To obtain the 
strange quark contribution to the singlet coupling, the 'loop' quark mass has been fixed 
to the strange quark massQ, while extrapolating in the light (valence) quark mass of the 
proton. 

We emphasize that the disconnected parts are of opposite sign to the connected ones, 
thus lowering the value of G\. Obviously, they agree within errors over the entire mass 
range. This indicates flavor symmetry of the disconnected amplitudes. The mass depen- 
dence of all the amplitudes is weak, and a linear extrapolation seems well justified. 



At the chiral limit, the authors of ref. [p.23|| obtain 

G\ = /\u- M = 0.638(54) + 0.347(46) = 0.985(25) (114) 

and 

G^ = Am + Arf + As = 0.638(54) - 0.347(46) - 0.109(30) = 0.18(10) , (115) 

where Ag = (P|g7375g|P). Note that G\ is significantly lower than the OZI estimate 
G\ ^ G\ = 0.579(25), even within a statistical uncertainty of about 50%. 

Unfortunately, the systematical errors are not well under control. As mentioned above, 
the comparison of the lattice estimate of G\ with the experimental result can be taken to 
give at least a rough impression on the size of systematical effects. From the inspection of 
eq. (|114|) and eq. (0) one estimates that those are of the order of 30%. 



■^^The lattice value of the strange quark mass rus is determined by the condition mKinis, mi)/mp(rni) 
0.64, where mi is the light quark mass. 
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Figure 11: Connected u-quark (a), connected d-quark (b) and disconnected contributions 
G\ as found by ref. I\12^] at a quark mass corresponding to mT,/mp ^ 0.745. The discon- 
nected part has been estimated with the (modified) volume source technique. The solid lines 
indicate the results of the fits according to eq. i\8^). 



A similar behavior of connected and disconnected amplitudes has been obtained by 
the authors of ref. MM. They find estimates G\ = 0.25(12), G\ = 1.20(10), and 
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Figure 12: Ref. ^12^J : Connected and disconnected amplitudes Aq = {Plqys'ynqlP) as a 
function of the bare quark mass in lattice units. The data has been multiplied by the 1st 
order tadpole improved renormalization constant Z\ = Z^'^ . 



G\ = 0.61(13), which are much closer to the experimental results. Since they have worked 
at a larger cutoff, a~^ ~ 2GeV compared to a~^ ~ l.SGeV in ref. ||123|| , the systematical 
errors might be indeed reduced. However, as mentioned above, the statistical sample which 
underlies this calculation is rather small, i.e. 1/10 of the sample used in ref, 
some caution in interpreting these findings is in order. 



1231. Thus, 



Nevertheless, the results found in quenched lattice calculations of G\ are encouraging, 
as they yield quantitative support for the assumption that the small value of G\ can be 
explained by negative contributions from quantum fluctuations. 

It will be interesting to see whether this remains valid in full QCD. 



3.4.2 Full QCD 

Recently, the SESAM collaboration has performed a full QCD lattice calculation of G\ 
||125|| . The analysis is based on 4 sets of 200 statistically independent gauge configurations, 
which have been generated in the Wilson discretization scheme with rtj = 2 dynamical 
fermions at /3 = 5.6. The cutoff a~^ = 2.3GeV corresponds approximately to that of a 
quenched calculation at /3 = 6.0 |]7D| . 

SESAM has applied the summation method, see eq. (PU]), and the standard insertion 
technique, c.f. section p72| , to calculate the connected contributions. We illustrate the 
quality of the corresponding signals in fig. O. Obviously the expected linear rise of 
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Figure 13: Full QCD, SESAM collaboration \12d^] : Connected contributions to G\ at 
different quark masses. The slope of the 'u quark' signals corresponds to the (spin averaged) 
amplitude s^(P|'U7mM75M|P), that of the 'd quark' signals to s^{P\d'yrnU'y5d\P) . The fits 
(range and value) are indicated by dashed lines. 



tdSUM 



it) is manifest for all four quark masses. 

The disconnected amplitudes have been determined with the PAM method, eq. ([91|). 
To calculate the numerator of R^^'^ SESAM has used the spin explicit stochastic estimator 



method, see section p.2| , with complex Z2 noise and 100 estimates per spin component. This 
procedure has been proven to be much more advantageous than the standard summation 
method in calculations of the disconnected contributions to ct^at [^ 
axial vector insertions are shown in fig. 



The signals for the 



TJ. Although a linear behavior can be identified, 
the statistical quality of the data is low. This is similar to the situation found in the 
quenched calculations with (un-improved) volume source and summation methods, see 
fig ~ 



11 



The apparent non-improvement of PAM and spin explicit stochastic estimator 
technique might be due to larger statistical fluctuations in full QCD compared to quenched 
QCD. As mentioned above, such a behavior has been recognized in the analysis of the light 
spectrum, where, at equal sample size and lattice setup, the statistical errors in full QCD 
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Figure 14: Full QCD, SESAM collaboration WMl : The ratio RA^^^t) with Ato = At = 1 
for the disconnected amplitudes Dq = {P{nsea)\{(ll3l5(l){i^sea)\P{i^sea)) o,t four sea quark 
masses. The fits (range and value) are indicated by solid lines. 



exceed those of quenched QCD by more than a factor of 2 [^. 

From the analysis of R^^'^ and R^^^^ at several quark masses one extracts the mass 
dependence of connected and disconnected contributions and extrapolates to the light 
quark mass. The data and the extrapolations are compiled in fig. IE. Obviously, the 
full QCD results are not very different from the quenched estimates of ref. ||123|| , which 
are shown in fig. O. In both cases one finds a weak mass dependence of the amplitudes. 
The disconnected insertions are less than zero, and their flavor symmetry appears to be 
maintained. 

At the light quark mass, the full QCD estimates for the axial vector coupling constants 
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Figure 15: Full QCD, ref. j\12^ : Extrapolation of the unrenormalized lattice amplitudes 
to the light quark mass. The fits are indicated by dashed lines, the results of the fits by 
bursts. The amplitudes are defined as follows: Cu = Sfj,{P{mq)\{ify^'y5u){mq)\P{mq)) con 
(spin averaged); Cd = s^{P{mg)\{d'y^'j5d){mq)\P{mq))con (spin averaged); Dq = 
{Pimq)\{q-fs-fr^q){mq)\P{mq)) disci D^ = (P(mg)| (573755) (m^)|P(mg))disc ■ 



read 
and 



G\ = Zi [Cu + Cd + 2Dq + Ds] = 0.20(12) 



C 



7NS 



A 



zr [Cu - Cd] = 0.907(20) , G\ = Z^^ [Cu + Cd] 



0.484(18) 



where the 1st order tadpole improved renormahzation 
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has been apphedQ. As in the quenched case at [3 = 5.7, the estimate of the triplet 
coupling differs from the experimental result by about 30%. This points to the presence 

^^SESAM has used a]v75(^) = 0.215, k = ki = 0.15846 for the hght quark insertions {Cu,Cd,Dq), and 
K = K, = 0.15608 for D,. 
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of sizeable systematic uncertainties due to insufficient renormalization and finite cutoff 
effects. However, the major source of uncertainty in G\ is given by the substantial amount 
of statistical fluctuations. 

Within these uncertainties one finds no significant difference between quenched and full 
QCD results. This might indicate that the disconnected contributions to the flavor singlet 
axial coupling of the proton are dominated by the gluonic properties of the vacuum. 

3.5 Summary 

Within statistical and systematical uncertainties of about 50%, lattice simulations consol- 
idate the assumption that the influence of QCD vacuum properties leads to a small value 
of G\. With the direct method, the latter are reflected in the disconnected contributions 
to the flavor singlet axial vector amplitude of the proton. Since quenched and full QCD 
calculations yield similar results, one might suspect that this phenomenon is dominated 
by gluons. 

Use of the topological method would allow to investigate the influence of the axial 
anomaly on the value of G\ in more detail. In particular one could study explicitely it's 
dependence on the topological sectors of the QCD vacuum. However, as we have pointed 
out above, this approach suffers from serious numerical problems, and one has not suc- 
ceeded yet to extract reliable estimates for G\ by application of this technique. 

Several strategies to improve this situation come in mind. With the direct method, 
the accuracy limiting factors are the statistical noise of the disconnected amplitudes and 
the neglect of higher order or non-perturbative contributions to Z^. The former can be 
reduced by use of larger lattices, where self averaging effects will help, higher statistics and 
adapted estimator techniques. To improve on the latter one needs, at least, a 2nd order 
calculation of Z^ in the framework of lattice perturbation theory. Even then, as can be 
seen from the findings of the QCDSF collaboration, one might have to fight against cutoff 
effects. Thus, a scaling analysis of G\ appears to be unavoidable. 

To apply the topological method successfully, one needs to generate (physically) large 
vacuum configurations at small cutoff. This would enable for a calculation of the correlation 
of the proton with the topological charge density at small quark masses and momentum 
transfer. Fortunately, a non-perturbative renormalization procedure for the topological 
charge is available. 

With current quenched or (n/ = 2) full QCD simulations, the strange quark contri- 
bution to G\ might not be estimated correctly . In a realistic setup one would have to 
include a strange sea quark. From the numerical point of view a lattice simulation with 2 
(degenerate) light quarks and 2 (degenerate) strange quarks appears feasible. 



52 S. Gijsken, Flavor singlet phenomena in lattice QCD 

4 The Mass of the r]' Meson and the U{1) Problem 

4.1 The Problem 

In a classical field theory, symmetries of the Lagrangian must be reflected by symmetries 
of the physical phenomena which such a theory is supposed to describe. In a quantum field 
theory, this is not necessarily true, since quantum fluctuations of the vacuum can spoil 
the classical symmetries. However, one then expects to obtain remnants of this symmetry 
breaking in the physical observables, e.g. the occurrence of (massless) Goldstone bosons in 
case of a spontaneously broken global symmetry and of massive gauge bosons and Higgs 
particles if a (local) gauge symmetry is violated. 

The QCD Lagrangian exhibits a f/(3) x f/(3) = SU{3) x SUa{3) x f/(l) x Ua{1) global 



flavor symmetry | 126 | in the chiral limit, niu = rud = rUs = 0. Most of the implications of 
this large symmetry on the physics of strong interaction has been understood already in 
the early days of current algebra and the quark model. 

The S'f/(3) flavor symmetry is reflected in the multiplet structure of the particle spec- 
trum[ [|127|| , whereas the S'f/A(3) axial flavor symmetry is supposed to be broken sponta- 
neously. The corresponding Goldstone bosons can be identified with the 8 light pseudo 
scalar mesons (7r°'^,i^°,Ar°,i^^ ,■?]), and the fact that those have non vanishing and non 
degenerate masses can be explained by an additional explicit breaking of the SU{?>) and 
S't/^(3) symmetries due to m^ ^ rrid ^ niu ^ 0. 

The f/(l) flavor symmetry predicts baryon number conservation, which is well realized 
in nature. Likewise, one would expect an ^ axial baryon number conservation', i.e. a 
parity doubling of baryon states, from the ^/^(l) symmetry. Since the particle spectrum 
does not exhibit such a parity doubling, one would suppose that this symmetry is broken 
spontaneously. However, there is no light pseudo scalar meson which could be identified 
with the corresponding Goldstone boson. The rj' meson has the correct quantum numbers, 
but it is simply too heavy. This is the so called U{1) problem of QCD: The Ua{^) symmetry 
of the QCD Lagrangian seems to be broken by quantum fluctuations, but apparently there 
is no remnant of this symmetry breaking visible in the particle spectrum. 



This problem is nearly as old as QCD itself. After Glashow ||126|| had noted in 1968 the 
cogency of a pseudo scalar meson with a mass of the order of the pion mass, it was shown 
by Weinberg ||128|| , using conventional soft pion methods, that this mass should indeed be 
smaller than y/S x mj^. 

Since then a natural solution to the U{1) problem has been achieved[^ by the study 
of the connection between the axial anomaly, topological excitations of the vacuum and 
instantons P,[l29i. 

The physical picture which emerges from this work is (roughly) as follows. The conser- 
vation of the ^a(I) symmetry current A^ , eq. (|69D , is formally violated by the presence of 



the axial anomaly 2n f^trG ^^G^'^ ■, c.f. eq. (0). One might wonder however, how such a 
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For reviews on the U(\) problem see [130-132 
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contribution can make a physical effect at all, if one assumes random quantum fluctuations 
of the gauge fields. In this case, a possible contribution from the axial anomaly should just 
average away. 

A possible solution to this problem came from G.'t Hooft |ll]], who showed that non- 
trivial topological excitations of the vacuum, instantons, can produce a non vanishing 
contribution to physical observables. 

In this scenario, the Ua{^) symmetry breaking would not be spontaneous and thus one 
would not expect to find a corresponding pseudo scalar Goldstone boson[^. Remnants of 
this new kind of symmetry breaking would however still be observable in nature. For ex- 
ample the rj' meson would acquire it's large mass from instanton contributions. If it were 
possible to switch off the topological vacuum excitations, the rj' meson should become light 
and would behave like an 'ordinary ' Goldstone particle. 

Indeed, such a behavior has been obtained in the framework of a l/N^ expansion, N^ 
being the number of colors, by Witten |jl2| and Veneziano |13|. They showed that in 



the limit Nc -^ oo, where the axial anomaly is absent, the masses of rj' and vr become 
degenerate. In first order 1/Nc this degeneracy is lifted and, in an SU{3) flavor symmetric 
world, one obtains 



,2 2 2 2n/x 



mr,'-m^ = mQ = —— . (119) 

J IT 

Here Uf denotes the number of flavors and x is the topological susceptibility defined in 
eq. (^), at /c^ = 0. Taking into account an explicit SU{3) flavor symmetry breaking and 
neglecting the effects of a (small) t] — rj' mixing, one would expect 

ml ~ m^, - {Am\ + 3m^ + ml)/S - (867MeV)2 , x ^ (182MeV)^ . (120) 

In summary, there is a lot of (theoretical) evidence, that the topological properties of 
the QCD vacuum are indeed likely to solve the U{1) problem. The question remains how- 
ever, whether this qualitative evidence holds also quantitatively. To answer this question 
one would first have to show that the large rj' mass, m^/ = 958MeV, is really a property of 
QCD. Secondly, one would like to see that the topological excitations in full QCD are such 
that they yield an additional contribution to m^/ of size mo — 867MeV. Last not least one 
would have to study the instanton content of QCD in detail to learn about the mechanism 
of how the rj' acquires it's mass. 

Being a non-perturbative method, lattice gauge theory is perfectly suited to study 
quark- ant iquark bound states like octet and singlet mesons. It has been argued in ref 's. 



'^''it has been argued frequently, see e.g. ref. [133|, that, even in the presence of the axial anomaly, it is 
still possible to define a conserved flavor singlet axial vector current. Thus, the U{1) symmetry would be 
broken spontaneously and there should still exist a corresponding Goldstone boson. Since this conserved 



current is not gauge invariant, the Goldstone boson could be realized as a Kogut-Susskind ghost |134| in the 
unphysical, gauge variant sector of QCD. The physical interpretation of such an unphysical construction 
however is by far not clear and requires at least additional assumptions about the nature of the vacuum 
pi. 
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131| , p.32|| that the contributions of topological vacuum excitations could be overlooked by 



standard perturbative expansions. Clearly, with lattice QCD this flaw is absent. Within the 
course of a Monte Carlo process, QCD vacuum configurations are explicitly generated^. 
Thus, one has direct access to topological properties of each single configuration. This 
enables, in principle, for a detailed study of the connection between instantons, topological 
excitations and the rj' mass. 

4.2 Lattice Technique 

The mass of a (full) QCD bound state can be extracted from the asymptotic behavior of 
the momentum zero projection of it's propagator. 

Let us consider, for example, a one flavor pseudo scalar meson with an interpolating 
operator P = qj^q. The time dependence of it's (momentum zero) propagator is given 
byS 

Gp(0,t)=^(pt(o,0)P(t,f)) = 5:i(0|pt|^)|V™"* '^ |(0|pt|no)|V"'"o*, (121) 
where uq is the lowest eigenstate of the QCD Hamiltonian with quantum numbers J^'" = 

0-+. 

Alternatively one can express Gp{0,t) in terms of quark operators 

Gp{0,t) = E((975g)(0,0)(g75g)(5^,t)) , (122) 

X 

or, performing all possible contractions of quark operators, in terms of quark propagators 
Ag, defined in eq.(^) 

Gp{0,t)= (123) 

5: {(Tr( A,At)(0,t)) - (Tr(A,75)(0,0)Tr(A,75)(t,t))} =C7p(0,t) -Dp(0,t) . 

X 

Thus, a one flavor meson propagator has always two contributions: The first term on the 
r.h.s. of eq.( |123| ), Cp, is called 'connected', since the meson states at time and time t 
are explicitly connected by quark lines. Correspondingly, the second term. Dp, is called 
'disconnected'. Note that this term violates the OZI rule. It is highly sensitive to the 
details of the vacuum structure since the correlation between the quark loops Tr(Ag75) at 
time and t is mediated by quantum effects (multiple gluons and internal quark loops). 
Both contributions are schematically depicted in fig.0. 
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One has to ensure of course that the Monte Carlo algorithm tunnels sufficiently often between different 



topological sectors |115] 

^^This is true for an infinitely extended lattice in time. For a finite lattice with time ex- 
tend T and (anti)periodic boundary conditions the asymptotic form is given by limt_,tx) G(0, t) = 
|(0|P'''|no)P [e"™"o* + e^™"o(^^*)]. For the sake of clarity we will however restrict the discussion to 
infinitely extended lattices. 
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The physical hght mesons are not made of one quark flavor. They are rather 2- flavor 
triplet (pions), 3-flavor octet (Kaons,ry) or 3-flavor singlet {r]') combinations of quarks. Just 
from combinatorics one expects the disconnected contributions to cancel each other largely 
for the non singlet combinations^, and to contribute additively to the singlet state. 

From this point of view, one can identify the connected contribution with the normal 
'Goldstone part' of the pseudo scalar mesons and the disconnected contribution with the 
anomalous 'topological part'. 



We could now proceed and calculate Gp according to eq. ( |123|) on a representative 



sample of gauge configurations. Then one would extract ?7i^/ from a fit to its asymptotic 



form, eq.( |121| ). However this might not be advantageous for several reasons. First, one 
would like to concentrate on the (anomalous) mass contribution from the disconnected 
part. Secondly, one should try to keep statistical errors, which arise mostly from the large 
fluctuations of he disconnected part, as small as possible. For both these reasons one favors 

to consider the ratio 

njDpiO,t) 

^^^^'^ - Cp{0,t) ' (^2^) 

rather than to Gp itself. The factor Uf (number of flavors) has been inserted in the 
numerator because a nj-flavor singlet meson has Uf contributions Gp, but n'j contributions 
Dp. 

In full QCD, the asymptotic time dependence of Gp and Gp is given by eq.( p,21| ). Thus, 
it follows 

R^c{t) = 1 - CeK'-^s)* (125) 

at large t. Here we have identified mass of the full singlet propagator (Gp) with m^' and 
the mass of the 'Goldstone part' (Gp) with rrig. Note that Rdc is sensitive to the difference 
of these masses. The 'Goldstone mass' mg vanishes in the chiral limit. Therefore, one can 
extract rriQ, see eq. ( |119| ), from an extrapolation to this point. 



The situation is slightly more complicated in quenched QCD. In this case eq.( |121| ) 
does not describe the asymptotic time dependence of Gp correctly since higher quark loop 
contributions to the singlet propagator are not taken into account. In the 1/Nc expansion 
of Witten and Veneziano [12,1^, which in first order leads to eq. ( p.l9| ) , the t]' propagator in 



momentum space, Gp, is described by an infinite sum over virtual quark loops, connected 
by an effective gluonic coupling ml 

p"^ + ml (p^ + mg)2 (p2 + mg)^ p^ + mg + mg 



This is shown in fig.|16[ The quenched approximation neglects virtual quark loops, which 
corresponds to a truncation of the geometric series in eq. ( |126| ) after the second term. 
Obviously, the first term can be identified with the connected contribution to the truncated 
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Up to a presumably negligible flavor symmetry breaking, they cancel exactly. 
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Figure 16: Contributions to the flavor singlet propagator in the scheme of Witten and 
Veneziano. 

propagator, whereas the second term belongs to it's disconnected part. A Fourier transform 
of both terms yields^ 

rrii , ml 



{RDc{t)r''' = A+-^t,A = -^. (127) 

2ms 2mi 



'■8 



Thus, under the assumption that the Witten- Veneziano expansion is vahd, one can de- 
termine rriri' even in the quenched approximation from a hnear fit to [Rncit))'^^'^"'- It is 
however by no means obvious, that quenched and full QCD calculations should yield the 
same value for m^/. Full QCD simulations with reliable statistics are therefore urgently 
needed. 

A peculiar situation arises in the partially quenched approximation, i.e. in the case 
where not all valence quarks have their counterparts in the quark sea. Such a setting can 
be useful if one wants to exploit gauge configurations with light sea quarks with respect 
to the strange quark sector of QCD or if one works with rif = 2 dynamical fermions in 
the Kogut-Susskind discretization. The form of the corresponding rj' propagator, i.e. the 
analogon to eq.( |126|) , has been worked out by Bernard and Golterman in the framework 
of the partially quenched chiral perturbation theory 



The numerical techniques which are used to calculate the connected and disconnected 
contributions to the rj' propagator, eq. (|123| ), are quite similar to those already explained 



in the previous sections. The connected part can be composed from solutions to eq.(|46|), 
and for the disconnected contribution one can use the estimator techniques discussed in 



sections |2.2| and |3.2| . We emphasize again that the computationally expensive and statisti- 
cally noisy part of this calculation is the disconnected one. In order to gain statistics one 
exploits the translational invariance of the propagator, G{0,t) = G{tQ,to + t), and performs 
the average over all timeslices Nt, Dp{0,t) = l/NtJ2to -Dp(^o,^o + t). 

We mentioned at the beginning of this section that the calculation of m^/ and the proof 
that it does not vanish in the chiral limit are not sufficient to proof that the 't Hooft scenario 
of the f/(l) problem is realized in QCD. One rather has to demonstrate that non-trivial 
topological excitations of the vacuum really exist in QCD, and that those excitations are 
responsible for the large mass of the r]' meson. 



^^Using the Witten- Veneziano expansion one can also estimate the amphtude C in eq.(125). One finds 
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One method to study the connection between topology and m^/ is to determine the 
topological charge Q, see eq.(|83|), in addition to m,,/ and mo, on a given set of gauge 
configurations^^ With this information in hands one can analyze the dependence of rrirj' 
(or equivalently of mo) on Q. A finding m,,/ = mg, (mo = 0), at Q = and m,,/ ^ rris, 
(mo 7^ 0), at Q 7^ would favor strongly the 'topological solution' of the f/(l) problem. 

Additional information on this connection can be extracted from the lowest eigenmodes 



of the Dirac operator. The Atiyah-Singer index theorem |136 



Q = n+ -n- (128) 

states that the topological charge Q in the continuum is given by the difference between the 
number of positive (ra+) and negative (n~) chirality eigenmodes \E'^ of the Dirac operator 
with zero eigenvalue 



fvi/± = ^M(9^_,^^A«)vl/± = zAvl/± , A = (129) 

75^+ = +^+ , 75^" = -^" . (130) 



Thus , if mo is determined by topology, one would expect that the related correlator, i.e. 
the disconnected part of eq.( |123| ), could be calculated using only the zero modes of the 
Dirac operator, whereas the latter should be insufficient for a complete determination of 
the connected part. 

Unfortunately, the Atiyah-Singer theorem does not hold (exactly) in a discrete formula- 
tion (Wilson or Kogut-Susskind) of QCD on a finite lattice. This is due to the fact that the 
fundamental property 75^ = — -P75 of the (anti hermitian) Dirac operator, is not fulfilled 
in the discretized case. It follows that, in principle, all eigenmodes of the lattice Dirac 
operator can contribute to the topological charge. 

It has been shown however, by Smit and Vink ||112|| and Itoh, Iwasaki and Yoshie ||137|| , 



that the Atiyah-Singer theorem holds still approximately on the lattice. Instead of taking 
into account only the zero eigenmodes one now needs the 'lowest' eigenmodes to make up for 
the topological charge. Since it is not a priori known how many of these 'lowest' eigenmodes 
are really needed, one has to study the infiuence of the number of lowest eigenmodes on 
the connected and disconnected contributions to the rj' propagator numerically. 

We mention that the determination of the spectrum of the Dirac operator on a (large) 
lattice represents a non-trivial numerical problem. Algorithms suited for this task are 



explained in ref's. [138-141 



4.3 Lattice Results 

The 'topological' scenario to explain the U(l) puzzle has been proposed more than 20 
years ago, but it is only in these days that the lattice method may become powerful 
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Lattice definitions and numerical methods to determine Q have been discussed in section |3.3| 



58 



S. Gijsken, Flavor singlet phenomena in lattice QCD 



enough to prove that the proposed solution is actually realized in QCD. This is mostly due 
to the substantial computational effort which is necessary to calculate the disconnected 
contributions to the rj' propagator with reliable statistics. 

Nevertheless there is already a handful of pioneering simulations in quenched QCD 
and, to an even lesser extend, in full QCD, which find promising evidence that the 't Hooft 
mechanism explains the large mass of the t]' meson. 



After early exploratory studies ||137| , |142| from the time where estimator techniques 
for the disconnected part of the rj' propagator were not yet available, the first quenched 
simulation with reliable statistics (200-300 gauge configurations) was performed by the 
authors of ref . ||143|| , using Wilson fermionsf^ on a 12^ x 20 lattice and a somewhat large 
lattice spacing, a ~ 0.14fm (a~^ = 1.45(3)GeV, P = 5.7). They employed the volume 
source technique, see eq's. (B3]),(B3|), to calculate the disconnected contribution. 



quenched 
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Figure 17: Left: Connected (Cp) and disconnected (Dp) contributions to the rj' propagator 
at a quark mass corresponding to niT^/mp = 0.59. Right: The ratio Rdc o,^ the same quark 
mass. The solid line indicates the result of a linear fit to Rdc ^^ the range 4 < t < 8. The 
data has been extracted from ref. ^J^. 



We illustrate the statistical quality of the corresponding signals in fig. 17, where we 
display connected and disconnected contributions to the rj' propagator as well as the ratio 
Rdc- The disconnected part is much more noisy than the connected one, as expected. 
The ratio Rdc is reasonably consistent with the predicted linear behavior in the quenched 
approximation, see eq. (|127|) . From this data, one can extract mo at a given quark mass. 



40 n 



They chose quark masses correspondmg to mTj/mp = 0.74, 0.70, 0.59, 0.42. Unfortunately, the results 
for the lightest quark mass were achieved with only 40 gauge configurations. 
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To mimic the physical situation of n/ = 3 quarks, one multiphes mo by -\/3. Clearly, such 
a procedure is justified only in a fiavor symmetric situation. 

The physical value for ruQ in this fiavor symmetric world, i.e. it's value at the light 
quark massQ can be read off from the result of a fit to mQ{mq). 



'0.55 



0.45 
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Figure 18: Ref. /\14^ -' Dependence of rriQ on the hare quark mass rriq. The solid line shows 
a linear fit to mQ{mq). The cross indicates the result of the fit at the chiral point. It's value 
in lattice units is rriQ = 0.518(25), which corresponds to m^ = 751(39)MeF. Note that all 
mo data have been multiplied by ^Jnj = y/S. 



As shown in fig. |T8|, the authors of ref. ||143|| find a linear increase of mQ{mq) with 



decreasing m„. Note that this is in contrast to the behavior of non-singlet pseudo scalar 



masses, where one obtains a linear decrease of m\, 



m„ 



At the chiral point. 



.m„ 



0), 



they get m-o = 751(39)MeV, where the number in brackets denotes the statistical error. 

One should, for various reasons, not be worried too much about the mismatch (2 — 
3cr) of this result with the physical estimate, see eq. (|120|) . As we mentioned above, the 



*^The light quark mass mi (in lattice units) can be determined by the condition mT^/mp = 0.1766. 
Alternatively one can extrapolate to the chiral point m; = 0. The difference between these methods with 
respect to toq is negligible. 
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disconnected part has been obtained with the volume source technique, which does not 
allow to monitor the cancellation of the unwanted gauge variant parts. Furthermore, the 
simulation has been performed at a coarse lattice spacing, but the influence of finite cutoff 
effects has not been taken into account. Moreover, one has assumed flavor symmetry with 
respect to mo, which might not be realized in nature. Last not least, the result has been 
achieved in the quenched approximation. Since vacuum contributions are essential in the 
determination of ttiq, this could introduce a systematic bias. For all these reasons one 
should expect that the true error on mo is actually much larger. 

The good news from the findings of ref. [|143|] is of course that mo is in the right range. 



Thus, there is a good chance that, once all the uncertainties listed above have been re- 
moved, one will indeed obtain the experimental 77' mass in full QCD lattice calculations. 

Given that, one can proceed to study the connection between vacuum topology and mg. 
The authors of ref. ||144|| have calculated the topological charge Q according to the standard 
definition^, see eq.( p.05| ), on the same set of gauge configurations as used in ref. ||143|| , at a 



quark mass corresponding to m.,^/mp = 0.59. Thus one can analyze the dependence of mo 



on Q for this setup. We display the result of ref. ||144|| in fig.|T9|. It appears that mo correlates 
to the value of Q: mo increases as Q increases. This indicates that the proposed solution 
of the U{1) problem is realized in QCD. Unfortunately, the authors of ref. |p.44|| have used 



an 'integerized' value of Q instead of a renormalized one, i.e. they set the lattice value of 
Q to the nearest integer instead of applying the renormalization procedure described in 



section 3.3. However, it is unlikely that this inaccuracy alters the conclusions of their work. 



In order to consolidate this result one needs to perform lattice simulations which avoid 
the insufficiencies discussed above. As this is, after all, a question of compute power and 
of the capacity of stochastic algorithms, this goal can be achieved only in small steps. 

Recently, Venkataraman and Kilcup ||146|| analyzed full QCD gauge configurations with 



2 flavors of Kogut-Susskind fermions with respect to the r/' mass. The configurations have 
been generated by the QCDSP collaboration | |147| , |148| on 16'^ x 32 lattices, at 3 different 



values of the sea quark mass, corresponding to mjr/mp = 0.704, 0.644, and 0.572, and at a 
lattice cutoff a~^ ~ 2GeV. The disconnected part of the r}' propagator has been calculated 



using a stochastic estimator method, and 'Wuppertal smearing' [|149|| has been employed to 
improve the ground state projection of the operators. In addition, the authors of ref. [|146|| 
have performed a quenched simulation at equal lattice spacing and volume. This allows 
for a direct comparison of quenched and unquenched results. Thus, this analysis has 



installed quite a number of improvements to alleviate the deficiencies of ref's. ||143| , |144 



Unfortunately, the size of the samples of gauge configurations used here is rather limited 
(34-79 configurations for full QCD and 83 configurations for quenched QCD). 

In fig.|3 W6 display the ratio Rdc both for the quenched and full QCD simulations. 
Obviously the expected t dependence of Rdci c.f. eq's. ( |125| ) and ( |127| ), is reasonably well 
satisfied by the data. 



^The definition is applied to 'cooled' configurations. The cooling method is explained in ref. |145] 
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Figure 19: mo in lattice units as a function of the 'integerized' topological charge Q\ at 
a quark mass corresponding to niT^/mp = 0.59. The solid line indicates the result for the 
entire ensemble of gauge configurations. The data has been extracted from ref. j\144\l - 



An extrapolation^ of the fit results at the various quark mass values to the chiral limit 
yields m^/ = 780(187)MeV for nj = 2 full QCD, and m^, = 891(101)MeV in the quenched 
approximation. As in ref. ||143|| , these values are obtained by correcting for the case of 3 
quark flavors. Within large statistical errors both results, quenched and unquenched, are 
consistent with experiment. This is similar to the situation found in lattice calculations of 
the axial flavor singlet coupling G\. It might indicate, within the given level of accuracy, 
that the vacuum contribution to m^/ is determined largely by gluonic properties, rather 



The authors of ref. |146| assume a hnear functional behavior of mQ{mq), in contrast to the finding of 



ref. [143|, where ttiq (not mp) increases hnearly with with decreasing niq. Unfortunately, the systematic 



uncertainty arising from this difference has not been taken into account. 
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Figure 20: The ratio Rdc for quenched {Nj = 0) and full {Nf = 2) QCD, as found by 
ref. ^14(^1 - The solid lines correspond to fits according to eq. l \123l ) and eq. ^12% ). 



than by fermion loops. 

We mention that this interpretation is supported by the findings of lattice calculations 
using the 'bermion' approach to QCD ||150|| . In this method one simulates QCD at negative 



numbers of flavor, nr < 0, where fermions are replaced by 'bermions', i.e. bosons with 



a fermion action, and then extrapolates the results to nj > 0. The authors of ref. ||151 



have applied this approach to the calculation of ttiq. As a result they find no significant 
dependence of mo on Uf in the range —8 < rif < 0. 



Finally, we comment on the numerical evidence for the connection between topological 
vacuum excitations and the rj' mass from a calculation of the low lying eigenmodes of the 
Dirac operator. According to the approximate validity of the Atiyah-Singer theorem on the 
lattice, the lowest eigenmodes refiect the topological properties of the vacuum. Thus, if the 
disconnected part of the rj' propagator in dominated by the latter, it should be possible to 
compose it exclusively from these contributions. The connected part, which represents the 
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Figure 21: Disconnected (left) and connected (right) contributions to the rj' propagator as 
found by ref. \14(\l in a full QCD (nj = 2) simulation with Kogut-Susskind fermions at a 
(lattice) quark mass ruqa = 0.01. The latter corresponds to mT,/mp = 0.572. 



'normal' Goldstone mode, should however not be obtainable from the lowest eigenmodes 
alone. 

The authors of ref. ||146[| have compared the connected and disconnected parts obtained 
from a full calculation with those constructed from the lowest 32 eigenmodes. We display 
their result in fig.|2T|. Obviously the (full) disconnected contribution is well represented by 
the eigenmodes, in contrast to the connected part. 



4.4 Outlook 

The final goal to confirm quantitatively that the 'topological' solution to the U{1) problem 
is realized in QCD has not been reached yet. There is however a lot of evidence that 
this solution is most likely. All results achieved so far with lattice methods point in this 
direction. 

From quenched and unquenched lattice simulations it appears, within the given level 
of uncertainty, that the mass of the r]' meson calculated in QCD is comparable to it's 
experimental value. The fact that both methods yield similar results may indicate that 
the topological contributions are dominated by the gluonic properties of QCD. 

On top of this, the connection to topological vacuum excitations has been established 
in quenched QCD by studying the dependence of mo on the topological charge Q, and in 
full QCD by a comparison between the r]' propagator calculated with all eigenmodes and 
the one composed from only the low lying states. 



To consolidate these findings one has to perform lattice calculations at high statistics 
and by application of (possibly improved) stochastic estimator techniques. This, as a first 
step, would enable for a reliable estimate of the statistical error. The next step is then 
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to vary lattice spacing and volume in order to perform a continuum extrapolation of the 
lattice results. If this is done within the framework of a 'conventional' full QCD simulation 
with 2 flavors and degenerate quark masses one could extract the QCD result for the r]' 
mass under the assumption of flavor symmetry with respect to rriQ and by neglecting rj — r]' 
mixing effects. 

In view of the fact that the results discussed above have been achieved on computers 
with a speed of less or about lOGfiops this aim is certainly within the range of upcoming 
TERAfiops machines. To include the effects of flavor symmetry breaking and to calculate 
reliably the influence oirj — rj' mixing effects one would have to repeat the whole procedure 
with nj > 3 dynamical quark flavors and non-degenerate masses. Without significantly 
improved stochastic algorithms, such a program is currently beyond reach. It should 
however be feasible to perform such a calculation with nj = 2 light and n/ = 2 heavy 
flavors at a fixed lattice spacing. This would help to estimate the size of these effects. 

Clearly, high statistics samples of QCD vacuum configurations will allow for a more 
detailed study of the connection between topology and m^/. Especially the role of instan- 
tond^, which we did not discuss here in detail, could be clarified. 



*For a recent review on instantons on the lattice see |152] 
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5 Final Remarks 

We said in the introduction that a major goal in studying flavor singlet phenomena is to 
learn about the properties of the QCD vacuum. So, what can we learn from the results of 
lattice QCD achieved so far in this context? The answer to this question embraces several 
aspects. 

First of all, we have seen that disconnected (vacuum) amplitudes indeed yield substan- 
tial contributions to flavor singlet processes. It turns out that they are suited in sign and 
order of magnitude to explain 'anomalous' phenomena like the small value of G\, the large 
mass of the rj' meson and the high value, compared to the octet estimate (Tq, of cr^j^. 

Secondly, the connection between topological excitations of the vacuum and the rj' mass 
could be established. This is evidence, that the f/(l) problem can be solved by the 't Hooft 
mechanism. Similarly, one flnds that the influence of the axial anomaly, which is (mostly) 
reflected in the size of the disconnected axial vector insertion of the proton, lowers the 
value of G\. 

Thirdly, comparing quenched and full QCD lattice results, one can estimate the effect 
of sea quarks on flavor singlet amplitudes. They seem to have some influence on the 
scalar insertions, e.g. on a^^N and y, as they lower the ratio of disconnected to connected 
amplitudes. Axial vector and pseudo scalar insertions, which are connected to the axial 
anomaly, do not seem to depend largely on the presence of sea quarks. 

Fourthly, flavor symmetry breaking effects in the disconnected (vacuum) contributions 
appear to be small. This would explain the success of lattice determinations of flavor non- 
singlet quantities like G\. 

Unfortunately most of the flndings mentioned here are still semi- quantitative results. 
At the current level of statistical and systematical reliability, the estimates of flavor singlet 
quantities from lattice calculations cannot be used yet for a precise quantitative comparison 
between experiment and theoretical QCD predictions. 

Clearly, one would like to have an idea of what it would take to reduce the uncertainties 
of the lattice results to the level of about 20%, which corresponds the accuracy currently 
achieved in experimental measurements of G\ and (TttAt- 

The various sources of uncertainty have been discussed in some detail in this review. 
It can be seen from this discussion that the error is dominated by the large statistical 
fluctuations of the disconnected amplitudes. These amount to 50 to 60%. To reduce the 
fluctuations to the experimental level of 20% just by 'brute force' one would need a factor 
of 9 in computer speed. 

The second largest error is due to cutoff effects and, closely connected to this, uncer- 
tainties in the determination of renormalization constants. From the lattice results of G\ 
one estimates a 20 to 30% contribution to the error. Thus, a rough flnite cutoff analysis 
would be sufficient. Current lattice calculations of flavor singlet quantities have been per- 
formed at a cutoff of a~^ ~ 2GeV. An increase of a~^ by a factor of 1.5 (a~^ ~ 2GeV) 
would require an increase of the number of lattice points by a factor of 1.5^ ~ 5. On top of 
this one has to take into account that it becomes more expensive to generate statistically 
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independent gauge configurations when going closer to the continuum hniit. For full QCD 
simulations, where the generation of gauge configurations consumes the largest fraction 
of computer time, the dependence of the auto correlation length on the cutoff is not well 
known. An additional factor of 2 in computer time should be sufficient however for a 
moderate increase of a~^. 

Thus, one estimates that, without any further algorithmic improvement, a factor of 
9 X 5 X 2 = 90 in computer power would be sufficient to reduce the uncertainties on flavor 
singlet calculations to the 20% level. 

According to the fact that the currently most advanced lattice calculations in this field, 
i.e. full QCD with Uf = 2 dynamical Wilson fermions, have been performed on (APE) 
computers with a sustained speed of 10 GigaFlops, one would need machines with a 1 
TeraFlops performance to achieve this goal. 



Since computers with a speed of a few hundred Gigaflops are already available [ |148| , |153 
we would like to close this review with the optimistic view that a 20% uncertainty in 
full QCD flavor singlet calculations can be realized within the next 5 years. Quenched 
simulations with that accuracy could be done right now. 
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